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Abstract 



The gravitino is a promising supersymmetric dark matter candidate, even 
without strict it!-parity conservation. In fact, with some small i?-parity vio- 
lation, gravitinos are sufficiently long-lived to constitute the dark matter 
of the universe, while the resulting cosmological scenario is consistent with 
primordial nucleosynthesis and the high reheating temperature needed for 
thermal leptogenesis. Furthermore, in this scenario the gravitino is unstable 
and might thus be accessible by indirect detection via its decay products. 

We compute in this thesis the partial decay widths for the gravitino in models 
with bilinear i?-parity breaking. In addition, we determine the neutrino sig- 
nal from astrophysical gravitino dark matter decays. Finally, we discuss the 
feasibility of detecting these neutrino signals in present and future neutrino 
experiments, and conclude that it will be a challenging task. Albeit, if de- 
tected, this distinctive signal might bring considerable support to the scenario 
of decaying gravitino dark matter. 



Zusammenfassung 

Das Gravitino ist auch ohne strikte i?-Paritat-Erhaltung ein vielversprechen- 
der Kandidat fiir supersymmetrische dunkle Materie. Tatsachlich sind Gra- 
vitinos bei leichter i?-Paritat-Verletzung geniigend langlebig um die dunkle 
Materie im Universum darzustellen, wahrend das daraus resultierende kosmo- 
logische Szenario mit primordialer Nukleosynthese und der hohen Reheating- 
Temperatur, die fiir thcrmischc Lcptogcncsc bcnotigt wird, vcrcinbar ist. Dar- 
iiber hinaus ist das Gravitino in diesem Szenario instabil und konnte daher 
durch seine Zerfallsprodukte fiir eine indirekte Entdeckung zuganglich sein. 

In dieser Diplomarbeit berechnen wir die partiellen Zerfallsbreiten des Gravi- 
tinos in Modellen mit bilinearer i?-Paritat-Brechung. Zudem bestimmen wir 
das Neutrinosignal astrophysikalischer Zerfalle von dunkler Materie, die aus 
Gravitinos besteht. SchlieBlich diskutieren wir die Moglichkeit diese Neutri- 
nosignale in heutigen und zukiinftigen Neutrinoexperimenten zu beobachten 
und kommen zu dem Schluss, dass es eine sehr herausfordernde Aufgabe 
ist. Dessen ungeachtct konnte dieses charakteristische Signal eine deutliche 
Unterstiitzung fiir das Szenario zerfallender Gravitinos als dunkler Materie 
darstellen. 
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Preface 



The nature of the dark matter (DM) in the universe is still one of the unsolved 
mysteries in modern cosmology. Its existence is strong evidence for physics 
beyond the Standard Model, but until now the only indication for dark matter 
is based on its gravitational interaction. 

One of the most favored solutions to the dark matter problem is the 
introduction of supersymmetry with conserved i?-parity. In this framework, 
the lightest supersymmetric particle (LSP) is stable and therefore a natural 
DM candidate, if it is also neutral and only weakly interacting. The most 
extensively studied case is where the LSPs are Weakly Interacting Massive 
Particles (WIMPs) and the dominant contribution to the dark matter in 
the universe. One of the most interesting candidates is, in fact, the lightest 
neutralino, since it can be tested in accelerator, direct detection and indirect 
detection experiments in the near future [1]. 

On the other hand, it is also possible that dark matter only interacts grav- 
itationally, and supersymmetry offers also candidates of this type. In local 
supersymmetry (i.e. supergravity) the gravitino arises as a natural candidate 
for dark matter. It is the superpartner of the graviton and has been the first 
supersymmetric dark matter candidate proposed [2]. The gravitino is also 
one of the most elusive dark matter candidates, since, as part of the gravity 
multiplet, all its interactions are suppressed either by the Planck scale (for 
the spin-3/2 component) or by the supersymmetry breaking scale (for the 
goldstino component). Nonetheless, a high reheating temperature after the 
inflationary phase in the early universe leads via thermal gravitino production 
to a relic abundance that can account for the observed dark matter density 
in the universe and is compatible with baryogenesis via thermal leptogenesis. 

Typically, the extremely weak interactions of the gravitino cause severe 
trouble for the standard cosmological scenario, since they make either the 
gravitino — if it is not the LSP — or the next-to-lightest supersymmetric par- 
ticle (NLSP) so long-lived that it decays during or after Big Bang Nucle- 
osynthesis (BBN). Therefore, the successful prediction of the light element 
abundances can be spoiled (see for instance [3]). 
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However, it was recently proposed that these problems are easily circum- 
vented if a small breaking of i?-parity is introduced [4]. In that case, the 
neutralino is too short-lived to play the role of dark matter, whereas the 
gravitino — if it is the LSP — can still have a sufficiently long lifetime, even 
exceeding the age of the universe by many orders of magnitude. Since the 
NLSP decays into Standard Model particles via i?-parity violating couplings, 
the cosmology reduces to the standard, non-supersymmetric case well before 
the onset of BBN. 

In collider experiments, late decays of the NLSP can be a signature of this 
gravitino LSP scenario. Another advantage of models with i?-parity violation 
is the fact that gravitino dark matter is no more invisible in astrophysical 
experiments: Even if the strongly suppressed couplings rule out an observa- 
tion in direct dark matter detection experiments, the tremendous number 
of gravitinos in the Milky Way and nearby galaxies compensates for the 
suppressed decay rate, possibly making gravitino decay products visible in 
indirect detection experiments. 

In fact, it has been found in [5, 6, 7] that gravitino dark matter with a life- 
time of ~ 10^^ s and a mass of ~ 150 GeV can account simultaneously for the 
anomalies in the extragalactic gamma-ray flux, suggested by the Energetic 
Gamma Ray Experiment Telescope (EGRET) [8], and in the positron frac- 
tion, suggested by the High Energy Antimatter Telescope (HEAT) balloon 
experiment [9]. 

Motivated by this coincidence, we compute in this thesis the neutrino 
signal in the same scenario, both as a consistency check and to find out 
whether in this scenario an anomalous contribution to the neutrino fiux may 
be expected in present and future neutrino experiments. While such signal is 
certainly much more difficult to disentangle from the background, it would 
provide significant support to the scenario of decaying dark matter, possibly 
consisting of gravitinos that are unstable due to bilinear i?-parity violation. 

The thesis is organized as follows: In the first chapter we will briefiy review 
Big Bang cosmology and the astrophysical evidence for dark matter in the 
universe. In Chapter 2 we will give a short introduction to supersymmetry and 
supergravity, and discuss more thoroughly the gravitino field, its interactions 
and its consequences for cosmology. In Chapter 3 we will introduce models 
with bilinear i?-parity breaking, discuss in detail the decay channels of the 
LSP gravitino in this framework and present the resulting neutrino spectrum 
from gravitino decay. 

In Chapter 4 we will start with the computation of the neutrino flux from 
extragalactic sources and from our own galaxy, considering also the effect of 
neutrino oscillations on the signal expected at the Earth. Then we will discuss 



2 



Preface 



the different neutrino baclcgrounds in the relevant energy range, compare 
them to our signal and propose strategies to disentangle the signal from 
the background. Subsequently, we will examine the prospects for a detection 
of the signal from gravitino dark matter decay in Super-Kamiokande and 
future neutrino detectors. Finally, we will present constraints on the gravitino 
lifetime from the non-observation of a neutrino signal, and then conclude in 
the last chapter. 

In Appendix A we introduce units with h = c = k = 1 and present a list 
of the physical and astrophysical constants used in the course of the thesis. 
In Appendix B we introduce the used notation and conventions. There, we 
also list numerous formulae used in the course of the thesis and derive some 
relevant identities. In Appendix C we present the Feynman rules for gravitino 
interactions and the detailed computations of the gravitino decay widths. In 
Appendix D we briefly present an example of the PYTHIA program used to 
simulate the fragmentation of the gravitino decay products. 

The results from this work have been published in [10]. 
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Chapter 1 



Dark Matter and the Early 
Universe 



In this chapter we want to review shortly the basic cosmological concepts that 
are relevant for the problem of dark matter. More comprehensive reviews on 
this topic can be found for instance in [1, 11]. 

1.1 Big Bang Cosmology 

In this section we will first introduce the dynamics of our universe and then 
highlight several important stages of its thermal history. 

Dynamics of the Universe 

The geometry of space-time is determined by the energy content of the uni- 
verse. This is expressed in Einstein's equations from general relativity 



Here, -R^j^ and R are the Ricci tensor and Ricci scalar, respectively, while 
Qf^u is the space-time metric and T^^, the energy-momentum tensor. To solve 
this set of coupled equations we have to assume symmetries of the universe. 
Measurements of the Cosmic Microwave Background (CMB) show that the 
universe is highly isotropic. In addition, galaxy surveys indicate that the 
universe is also homogeneous on large scales ((9(100) Mpc). 

The most general space-time metric compatible with isotropy and homo- 
geneity is the Friedmann-Robertson-Walker metric. The line element reads 





g^udx^dx" = ds^ = dt^ - a^{t) 



dr"^ 





1 — fcr^ 
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where a{t) is the scale factor, r, 9 and are the comoving spatial coordinates 
and the constant k characterizes the spatial curvature of the universe: k — —1 
corresponds to an open, A: = to a fiat and /c = +1 to a closed universe. 
Substituting dx = dr/ \J\ — kr^, the line element can be written as 

ds^ = df - a\t) [dx^ + S^{x) {dO^ + sin' 6 , (1.3) 

where S{x) — sinhx for k — —1, S{x) = X for /c = and S{x) — sinx for 
k = +1. 

A simplifying assumption — compatible with the above discussed symme- 
tries of the universe — is that the matter and energy content of the universe 
can be described by a perfect fluid. The energy-momentum tensor for a per- 
fect fluid in its rest frame reads 



/+g \ 





\ -pj 



(1.4) 



Solving Einstein's equations with the above assumptions results in the 
Priedmann equation 



= H ^—^^g,- (1.5) 



a / 6 — ' a 

and the acceleration equation 

^ = i^ + i^^ = -^E(^^ + 3A)- (1-6) 

Cl o 

i 

These equations determine the dynamics of the universe. Here we introduced 
the Hubble parameter H = a/a that characterizes the expansion rate of the 
universe. Prom the Priedmann equation and the acceleration equation we can 
derive the continuity equation 

A + 3i^(^i+pO -0, (1.7) 

which is equivalent to the covariant conservation of the energy-momentum 
tensor. 

There are several kinds of energy content in the universe, characterized 
by the equation of state 

Pi = WiQi. (1.8) 
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Radiation and relativistic particles have = 1/3 whereas pressureless mat- 
ter or dust (i.e. non-relativistic particles) has Wm = 0. A cosmological con- 
stant A can be described by an energy component with negative pressure 
wa = -1. 

The energy of photons and other relativistic particles decreases during 
their propagation through an expanding universe. This is expressed in terms 
of the redshift parameter z that is defined to be 

^ = ^-1 + -, (1.9) 

where Ae and tte arc the wavelength and scale factor at emission, and Ao and 
Qq are the observed wavelength and the present-day scale factor, respectively. 

From the continuity cqiiation we can derive the dependence of the energy 
density on the redshift parameter 

g,{z) = g^l + zf^'^'^ . (1.10) 

The energy density of non-relativistic matter decreases with (1 + z)"^ due to 
the dilution of the number density with the expansion of the universe. By 
contrast, the energy density of relativistic matter decreases with an additional 
factor of (1 -I- 2;) because of the energy redshift in an expanding universe. The 
cosmological constant is equivalent to an intrinsic energy of the vacuum and 
is independent of the dynamics of the universe. 

We can rewrite the Friedmann and acceleration equations using the den- 
sity parameter 

a(.) = ^ (1.11) 

that gives the energy density with respect to the critical density gdz) = 
SH^/^SttGn), that corresponds to a spatially flat universe. The present-day 
critical density is given by 

OTj2 

= — ~ 1.05 X lO-^/i^GeVcm-^ (1.12) 

where the Hubble constant Hq is the present-day Hubble parameter and is 
usually parameterized as 

//o = 100/ikms-^Mpc-\ (1.13) 

with h ~ 0.7. With the above definitions the Friedmann equation and the 
acceleration equation can be rewritten as 

k k 

1 = ^tot{z) 5^ = ilr{z) + ^m{z) + ^a{z) - 



ilr (1 + Zy + + ZY + - ^j-^ (1 + Z) 



aim 
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-\ fitot(^) (1 + 3 w,s{z)) = -\Y. ^*(^) (1 + 3 w^i) • (1-15) 



Thus, the present-day Friedmann equation gives a cosmic sum rule 

k k 



Cosmological observations suggest that the universe is spatially flat. Addi- 
tionally, the present-day radiation density is negligible, giving the relation 



l~fitot^fi„^ + fiA. (1.17) 



Inflation 



Standard cosmology requires very specific initial conditions in order to ex- 
plain the observed cosmological parameters. The main two problems are the 
flatness problem and the homogeneity problem. 

Differentiating equation (1.14) with respect to the time, we obtain 

"^^^'^/^^ = -2 4- (fitot - 1) = (1 + 3^eff) HQ,,, (fitot - 1) . (1.18) 
dt art 

Since the Hubble parameter is always positive in an expanding universe, we 
see that Qtot departs from 1 in a matter or radiation dominated universe, if it 
is not exactly 1 in the beginning. Thus, in order to obtain the observed value 
fitot — 1 today, the initial value must be extremely fine-tuned. This dilemma 
is known as the fiatness problem. 

Cosmic microwave background observations indicate that the universe 
was highly isotropic before structure formation. However, the observed CMB 
sky is many orders of magnitude larger than the causal horizon at the time 
of photon decoupling. Thus, the homogeneity of the temperature could not 
be achieved by physical interactions. Instead, it could only be achieved by 
extremely fine-tuned initial conditions. This dilemma is known as the homo- 
geneity problem. 

Both problems can be solved by the introduction of an inflationary phase 
(where ~ — 1) that lasts for about 60 e-folds. This means that the scale 
factor grows by a factor of in this phase. During infiation, Qtot = 1 
becomes an attractor as we can see from equation (1.18). Therefore, the 
universe can arrive at Qtot — 1 regardless of the initial conditions and stay 
close to that value until today. The isotropy of the CMB sky can also be 
explained: The entire observed universe had initially been a small causally 
connected region and had expanded tremendously during infiation. 
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Such an inflationary phase can be reahzed by a scalar inflaton field that 
enters a so-called slow-roll phase. Apart from solving the above issues, infla- 
tion theories predict large-scale density perturbations that arise from quan- 
tum fluctuations of the inflaton field. These are observed in the form of 
temperature anisotropics in the CMB and finally lead to the formation of 
structures like galaxies and stars in the universe. 

After the inflationary phase, the density of all particles that initially were 
in the universe is diluted. However, the decay of the inflaton field at the 
end of the inflationary phase transfers its energy density into a hot thermal 
plasma of elementary particles. This process is known as the reheating of the 
universe, and the equilibrium temperature of the thermal plasma is therefore 
called the reheating temperature. After this phase the universe is described by 
standard thermal cosmology. Energies in the early universe are usually given 
according to the characteristic temperature of the thermal plasma. Due to 
the adiabatic expansion of the universe the plasma temperature decreases as 



according to the redshift of relativistic particles. Here, Tq is the present-day 
CMB radiation temperature. 

Baryogenesis via Thermal Leptogenesis 

A crucial question in cosmology and particle physics is why there is more 
matter than antimatter in the universe. This problem shows up in the baryon- 
to-photon ratio r/ = {rib — ni)/n^ = rib/n^ that is different from zero. In 
order to generate a baryon asymmetry it is necessary to satisfy the Sakharov 
conditions [12]: 

• Baryon number {B) violation, 

• C-symmetry and CP-symmetry violation, 

• Departure from thermal equilibrium. 

One of the models proposed to solve the problem of baryon asymmetry is 
baryogenesis via thermal leptogenesis [13]. In this model baryon asymmetry 
is generated from a lepton asymmetry. A nonvanishing lepton number L can 
be converted into a nonvanishing baryon number through non-perturbative 
sphaleron processes. These processes violate the linear combination B + L 
but conserve B — L.li the sphaleron processes are in thermal equilibrium in 
the reheating phase of the universe, they lead to the relation [14] 



T = To (1 + , 



(1.19) 



B 



2A + ANh 
66 + 13 Nh 



(B-L), 



(1.20) 
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where Nh is the number of Higgs doublets and a theory with three fermion 
generations is assumed. In the minimal supersymmetric Standard Model with 
two Higgs doublets (see Section 2.1) one has i? = 8/23 (5 — L) = —8/15 L. 

In thermal leptogenesis the needed lepton number L can be created in CP 
violating out-of-equilibrium decays of heavy right-handed Majorana neutri- 
nos. This mechanism is closely related to the problem of neutrino masses, 
since heavy right-handed Majorana neutrinos can also explain small non- 
vanishing masses for the light neutrinos via the seesaw mechanism. The ac- 
tual observation of nonvanishing neutrino masses in the last years strongly 
supports the existence of heavy right-handed neutrinos and therefore also the 
mechanism of thermal leptogenesis. 

To achieve the observed baryon asymmetry, the model of baryogenesis via 
thermal leptogenesis needs a high reheating temperature in the early universe 
oiTR>WGeV [15, 16]. 



Primordial Nucleosynthesis 

Big Bang nucleosynthesis takes place at T ~ 1-0.1 MeV and is thus based on 
well understood Standard Model physics. Moreover, BBN offers the deepest 
reliable probe of the early universe. The predictions of the abundances of the 
light elements D, '^He, ^He and ''Li are very sensitive to the physical conditions 
at that temperature. Thermal equilibrium of neutrons and protons is achieved 
through weak interactions like 

nue-^pe', ne^^pue and n^pe~Ue, (1-21) 



and leads to a neutron-to-proton ratio of 



rir. 



(1.22) 



where we neglected the chemical potential. Especially the number of relativis- 
tic particle species (e.g. the number of light neutrino species) and the baryon- 
to-photon ratio 7] = nb/n^ determine the freeze-out time of weak interactions 
and thereby fix the initial neutron-to-proton ratio to be nn/up ^1/6. This 
ratio slightly decreases to about 1/7 due to neutron decay until the neutrons 
are stabilized in bound states. Regardless of the detailed interaction pro- 
cesses, virtually all neutrons combine with protons to form ^He. The relative 
abundance by weight of ^He can then easily be estimated: 

Y = ^ ^ = !^L^ ~ 25 %. (1.23) 

Qp + Q^c np + rin 1 + n„/ Up 
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Figure 1.1: Q\ — plane of the concordance model of cosmology (left). 
The combination of several astrophysical observations suggests — 0.75 and 
~ 0.25. BBN predictions of the light element abundances (right). All bands 
show the 2 a uncertainties. The boxes show the primordial abundances derived 
from astrophysical observations (2 a statistical and 2 a statistical + systematic 
errors, respectively). Both figures are taken from [17]. 



The complete calculation of the light element abundances involves all the 
details of nuclear interactions and is able to predict the abundances at the 
10"^° level. 

The Standard Model particle content and a baryon-to-photon ratio r] ~ 
6 X 10~^°, as determined from the CMB measurements [17], yield abundances 
of the light elements that are in good agreement with data from astrophysical 
observations (see Figure 1.1). 

The agreement of the predictions and measurements of the light element 
abundances from BBN constrains deviations imposed by physics beyond the 
Standard Model. 

Cosmic Microwave Background 

The cosmic microwave background is a relic from the time when the photons 
decoupled from the thermal plasma of electrons and light elements at T ~ 
0.25 eV (i.e. z ~ 1100). After the temperature of the photons dropped below 
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the energy to ionize hydrogen, the universe became transparent for them. 
The Cosmic Background Explorer (COBE) satelhte mission found that the 
CMB is highly isotropic and corresponds to an almost perfect black body 
radiation spectrum with a temperature of Tq ~ 2.7 K (i.e. Tq ~ 2.3x10"^ eV). 
Furthermore, COBE observed temperature anisotropies of the CMB at the 
10"^ level. 

The Wilkinson Microwave Anisotropy Probe (WMAP) satellite mission 
investigated the CMB temperature anisotropies in detail. Expansion of the 
temperature anisotropy map into spherical harmonics 

AT °° I 

— {6, <P) = J2J2 (^in^YimiO, <P) (1-24) 

1=2 m=-l 

gives the CMB power spectrum l{l + l)C;/(27r) in terms of the multipole 
moment / with 

Ci = {\airnf) = X] \aimf . (1.25) 

m=—l 

Using a cosmological model with a limited number of parameters, it is possible 
to obtain best-fit values for the cosmological parameters from the observed 
CMB power spectrum. Generally, one uses a ACDM model that assumes a 
cosmological constant A and Cold Dark Matter (CDM, i.e. particles that 
were already non-relativistic before structure formation) as the dominant 
dark matter component. Some of these parameters are given in Appendix A 
(see also Figure 1.1). In particular, from the position of the first acoustic 
peak one finds that the universe is approximately flat {Qtot — 1)- 

1.2 Dark Matter 

In the following we will discuss the astrophysical evidence for dark matter in 
the universe, briefly touch upon the proposed particle physics candidates for 
dark matter and introduce two astrophysical particle dark matter detection 
techniques: Direct and indirect detection of dark matter. 

1.2.1 Evidence for Dark Matter 

Astrophysical observations suggest the existence of dark matter: There is ev- 
idence from observations on galactic scales up to cosmological scales through 
the gravitational force. On the other hand, there is no experimental proof of 
DM at microscopic scales yet. In this section we summarize the astrophysical 
DM evidence following the review in [1]. 
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r (kpc) 



Figure 1.2: Typical rotation curve of a galaxy (left). The solid line shows the 
fit to the data (dots with error bars) and the dashed, dotted and dash-dotted 
lines show the decomposition into the contributions from the luminous disk, 
gas and the dark matter halo, respectively. Figure taken from [18]. Shapes of 
the halo density profiles listed in the text (right). All profiles except for the 
isothermal are singular at the galactic center. 



Galactic Scale A very convicing evidence for dark matter on galactic 
scales comes from the observation of rotation curves of galaxies, i.e. the cir- 
cular velocity distribution of stars and g function of the distance to 
the galactic center. 

From Newtonian dynamics we expect the circular velocity to be given by 



^(''^ = y ' (1-26) 

with M(r) = 47r JJ" g{r)r'^dr. In the outer regions of a galaxy (where there 
is no luminous matter anymore) one would expect the velocity to fall off as 
V oc r~^/^. However, observations of galactic rotation curves show that the 
velocity remains constant even far beyond the luminous disk (see Figure 1.2). 
This can be explained by the existence of a spherical halo of dark matter with 
a density profile Qhaio oc 1/r^ in the outer regions. Different halo profiles can 
be parameterized in the following way: 

Qhalo{r) = ^ a.{f3-y)/a ' i^-'^'^) 

{r/rc) [I + [r/r^) 

Profiles obtained from N-body simulations include the profiles of Navarro, 
Frenk and White (NFW) [19], Moore et al. [20] and Kravtsov et al. [21]. Ad- 
ditionally, we list two profiles including adiabatic compression (i.e. including 
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Profile 


a 




7 


Tc (kpc) 


^0 (GeV cm ^) 


NFW 


1.0 


3.0 


1.0 


20.0 


0.26 


NFWc 


0.8 


2.7 


1.45 


20.0 


0.16 


Moore 


1.5 


3.0 


1.5 


28.0 


0.059 


Moorec 


0.8 


2.7 


1.65 


28.0 


0.064 


Kra 


2.0 


3.0 


0.4 


10.0 


0.57 


iso 


2.0 


2.0 





3.5 


2.1 



Table 1.1: Parameters for the halo models listed in the text. The normaliza- 
tion is given according to the local halo density in the Milky Way. 

the effect of the baryonic gas in the halo): The NFW compressed and the 
Moore compressed profile. For comparison we also list the simple isothermal 
profile. The parameters of these profiles are presented in Table 1.1 and taken 
from [1, 22]. The density has been normalized using a local dark matter den- 
sity of Qioc = 0.3GeVcm~^ for the Milky Way [1]. While the slope in the 
outer regions is strongly constrained by observations, the slope of the halo 
density profile in the inner parts of galaxies is not very well known: All halo 
profiles from N-body simulations have a singular behavior at the galactic 
center, whereas the isothermal profile is finite (cf. Figure 1.2). 

It has been proposed that part of the galactic DM halos could be com- 
posed of non-luminous Massive Compact Halo Objects (MACHOs). These 
could be for instance stellar-mass black holes, feint neutron stars, white 
dwarfs or stars, or Jupiter-sized planets. There have been searches for these 
objects through the microlensing effect, where several collaborations moni- 
tored the luminosity of millions of stars in the Large and Small Magellanic 
Clouds for several years. However, these experiments found that MACHOs 
can only contribute a subdominant part of the galactic DM. Thus, we need 
non-baryonic dark matter to explain the galactic dynamics. 

Scale of Galaxy Clusters Evidence for non-baryonic dark matter comes 
also from the scale of galaxy clusters. In order to indicate DM in clusters, one 
has to show that the cluster mass inferred from the luminous components is 
smaller than the mass obtained independently using the gravitational force. 
Several methods provide this possibility. 

Application of the virial theorem to the observed velocity dispersion of 
galaxies in clusters allows an estimation of the cluster mass. Observation of 
the X-ray spectrum of the hot gas inside clusters gives the temperature of the 
gas, which is correlated to the cluster mass under the assumption of the gas 
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being in hydrostatic equilibrium. Another method is gravitational lensing: 
Light from distant sources is distorted due to the mass of galaxy clusters 
in the line of sight. Analysis of a combination of many background sources 
allows to determine the shape of the potential well of the cluster. 

Comparison of the cluster mass estimated from luminous matter and gas 
to the mass determined using the above listed methods gives a clear evidence 
for DM in galaxy clusters. 

Cosmological Scale In contrast to the previously discussed observations, 
the measurement of the cosmic microwave background can be used to deter- 
mine the total amount of dark matter in the universe. The values found for 
the matter density and the baryon density of the universe in analyses of the 
CMB power spectrum strongly hint on the existence of a non-baryonic dark 
matter. 

The combination of all astrophysical observations gives striking evidence 
for a cold non-baryonic dark matter component of the universe. The value 
of Qm from cluster observations, CMB measurements and supernova data 
(concordance model, see Figure 1.1) together with the baryon density Qb 
from BEN predictions and CMB measurements leads to an energy density of 
non-baryonic dark matter of Qnivih"^ = ^mh"^ ~ VL^h'^ ~ 0.1 (cf. Appendix A). 

1.2.2 Particle Dark Matter Candidates 

We have seen that the dark matter in the universe must be non-baryonic. The 
observations also imply that the dark matter is electromagnetically neutral 
and has interactions comparable to the strength of weak interactions or even 
smaller. 

The Standard Model (SM) of particle physics offers only one candidate for 
particle dark matter: neutrinos. They are neutral and weakly interacting, and 
due to their small mass they would contribute to hot dark matter (i.e. they are 
still relativistic at the time of structure formation). Experimentally, however, 
simulations of structure formation favor cold dark matter in order to be in 
agreement with the observed large scale structure. In addition, it was found 
that the three light neutrinos cannot amount to a significant energy density: 
0.0005 < VL^K^ < 0.023 [17]. Thus, SM neutrinos do not play a dominant role 
for the dark matter problem. However, there is still the possibility that sterile 
neutrinos, that are neutral with respect to the SM gauge groups, account for 
the dark matter of the universe. 

Apart from that, several theories beyond the Standard Model propose 
natural dark matter candidates (see for instance the reviews in [1, 23]). The 
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most prominent and most extensively studied is the lightest neutralino of 
supersymmetry. With conserved i?-parity — if the neutralino is the lightest 
supersymmetric particle — it would be absolutely stable. Additionally, the 
neutralino can have naturally a relic density on the correct order of magni- 
tude. Other supersymmetric candidates are the superpartner of the neutrino, 
the sneutrino, and the superpartner of the axion, the axino. Sneutrino dark 
matter has already been ruled out by direct dark matter detection experi- 
ments, while the axino is still a viable dark matter candidate. 

In supergravity the gmvitino arises as a new dark matter candidate with 
extremely weak interactions. Thus, if it is not the LSP, it is very long-lived 
and thereby the origin of many problems in standard cosmology. On the other 
hand, if the gravitino is the LSP, it is stable due to i?-parity conservation 
and the NLSP becomes long-lived, causing similar cosmological problems. 
However, this work will deal with the consequences for indirect detection of 
gravitino dark matter in models, where i?-parity is slightly violated. In these 
models, all the cosmological obstacles can be evaded (see Section 2.3). 

Even though there are a lot more particle dark matter candidates apart 
from supersymmetry, we will restrict to this list since we are primarily inter- 
ested in the gravitino in this work. 



1.2.3 Direct and Indirect Detection of Dark Matter 

Besides collider experiments that can determine the mass spectrum of par- 
ticles beyond the Standard Model, there are two astrophysical techniques to 
detect dark matter particles: Direct and indirect detection of dark matter 
(see for instance [1]). 



Direct Detection 

Direct detection experiments try to detect particles from the dark matter 
halo that cross the Earth. They observe the recoil of target nuclei induced 
by elastic scatterings of weakly interacting massive particles off those nuclei. 
The expected signal is determined assuming that the dark matter particles 
are distributed in the halo according to an isothermal profile and have a 
Maxwell-Boltzmann velocity distribution with a characteristic velocity of 
Vq ~ 270kms-i [1]. 

Gravitinos have much weaker interactions than WIMPs and thus no signal 
from gravitino dark matter can be observed in direct detection experiments. 
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Figure 1.3: Comparison of the angular profile of the flux from decaying DM 
(solid line) and annihilating DM (dashed line) for a NFW halo profile. Here, 
6 is the angle to the galactic center. The fluxes are normalized to their values 
at the galactic poles {6 = ±90°). Figure taken from [24]. 

Indirect Detection 

The method of indirect dark matter detection uses the observation of cosmic 
rays produced in dark matter annihilations or decays. The cosmic ray flux 
originating from dark matter is proportional to the annihilation rate of the 
dark matter particle or — if it is unstable — to its decay rate. Annihilations 
require the collision of two dark matter particles, so the resulting flux de- 
pends on the square of the dark matter density. The flux from decays, on the 
other hand, is proportional to the density. Fluxes from annihilations therefore 
depend strongly on the dark matter distribution. 

A comparison of the fluxes J from annihilating and decaying dark matter 
in a galactic halo with NFW density distribution is shown in Figure 1.3. 
While the signals from the galactic center and anticenter differ by less than 
two orders of magnitude for the case of decaying dark matter, they differ by 
more than four orders of magnitude for annihilating dark matter. The exact 
numbers depend on the angular resolution used to regularize the flux from 
the galactic center, since the NFW profile is singular at that position. 

Cosmic rays from annihilating dark matter are therefore most likely ob- 
served from very dense regions like the galactic center. Other astrophysical 
objects like the Sun or the Earth could also amplify the dark matter density, 
since weakly interacting dark matter particles are captured in the center of 
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these objects via energy loss through scattering off nucleons. However, only 
neutrinos can escape these dense objects and could be observed. 

In the case of gravitino dark matter there is no accumulation inside as- 
trophysical objects since the gravitational interaction is too weak. Thus, we 
have a gravitino distribution according to the galaxy density profile. There 
is no chance to observe gravitino annihilations, since these are suppressed by 
higher orders of the Planck scale. So, if the gravitino is stable, there is no 
observable signal in indirect detection experiments. On the other hand, if the 
gravitino is unstable, it could be possible to detect signals from its decay. 
In this thesis wc want to discuss the prospects for the indirect detection of 
neutrino signals from the decays of an unstable gravitino. 
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Supersymmetry, Supergravity 
and the Gravitino 



In this chapter we want to give a brief introduction to the relevant top- 
ics of supersymmetry and supergravity and then discuss in more detail the 
gravitino. We will mainly follow the reviews on supersymmetry and the in- 
troductions to supergravity given in [17, 25, 26, 27, 28]. For the discussion of 
the gravitino field we refer to [29, 30, 31]. 

2.1 Supersymmetry 

Supersymmetry (SUSY) is a generalization of the space-time symmetries of 
quantum field theory (QFT) that relates bosons and fermions. It introduces 
new fermionic generators Q that transform fermions into bosons and vice 
versa: 



This is a nontrivial extension to the Poincare symmetry of ordinary QFT 
and its structure is highly constrained by the theorem of Haag, Lopuszanski 
and Sohnius [32]. 

If SUSY were an exact symmetry of nature, particles and their super- 
partners would be degenerate in mass. However, since no superpartners have 
been observed yet, SUSY must be a broken symmetry. 

Although not yet confirmed experimentally, there are several theoretical 
motivations for interest in this additional symmetry. The first is the hierarchy 
problem of the Standard Model. This problem stems from the huge difference 
of the electroweak scale (C(IOO) GeV) and the (reduced) Planck scale 



Q |boson) ^ |fermion) , 



Q |fermion) ~ | boson) . 



1 



2.4 X 10'^ GeV, 



(2.1) 



Mpi 
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Figure 2.1: Renormalization group evolution of the inverse gauge couplings 
in the Standard Model (dashed) and the MSSM (solid). The sparticle mass 
thresholds are varied between 250 GeV and ITeV. Figure taken from [25]. 

where gravitational interactions become comparable in magnitude to gauge 
interactions. The only scalar particle in the Standard Model, the Higgs boson, 
receives quadratic radiative corrections to its mass due to fermion loops. 
These quadratic divergences can be cancelled by the contributions of the 
bosonic superpartners to the radiative corrections. However, this cancellation 
works only for softly broken supersymmetry. If the Higgs mass is to be at the 
TeV scale, the soft SUSY breaking parameters can be no larger than a few 
TeV. 

The second motivation is the unification of gauge couplings = g^/Air, 
a = 1, 2, 3, where gi = ^^/5/3g' and g2 = g are the electroweak coupling con- 
stants and g^ = gs is the strong coupling constant of the unbroken Standard 
Model gauge group 

SU{3)c X SU{2)l X f/(l)y . 

The renormalization group equations, that determine the evolution of the 
gauge couplings, depend on the particle content of the theory. In the Standard 
Model the gauge couplings do not unify. By contrast, with the altered particle 
content of a supersymmetric theory at the TeV scale, the gauge couplings 
unify at the unification scale Mu ^ 2 x 10^^ GeV (see Figure 2.1). 

The third motivation is that SUSY provides a promising candidate for 
particle dark matter: The lightest supersymmetric particle. With conserved 



19 



Chapter 2. Supersymmetry, Supergravity and the Gravitino 



Name 


Gauge 


bosons 


Gau} 


2;inos 


(5f/(3),, SU{2)l)y 


B boson, bino 
W bosons, winos 
gluons, gluinos 


.(1) 

4 (2) a 


= B, 

= wi 


^(3)" 


= B 
= W" 

= r 


(1, 3)o 
(8, l)o 



Table 2.1: Gauge supermultiplets of the MSSM. 



i?-parity the LSP is absolutely stable. In order to constitute the dark matter 
of the universe, the LSP must be colorless and electromagnetically neutral. 
In that case, it only interacts gravitationally and — if it is not neutral with 
respect to weak interactions — through weak elastic scatterings with Standard 
Model particles. 

The Minimal Supersymmetric Standard Model 

The Minimal Supersymmetric extension of the Standard Model (MSSM) in- 
troduces superpartners of the Standard Model particles as well as two Higgs- 
doublets and their superpartners. Standard Model particles and their super- 
partners are combined into supermultiplets. The gauge supermultiplets (or 
vector supermultiplets) consist of a spin-1 vector boson A^, a spin-1/2 Ma- 
jorana fermion A"^ and a scalar auxiliary field D", where a labels the gauge 
group generators. In particular, we have the electroweak gauge bosons and 
the corresponding fermionic gauginos as well as the gluons and the fermionic 
gluinos (see Table 2.1). The table gives the notation of the particles in La- 
grangians and Feynman rules, and their transformation properties under the 
Standard Model gauge groups. 

Chiral supermultiplets consist of one complex scalar (f), a two-component 
chiral fermion x and an auxiliary scalar field F. Namely, there are three gen- 
erations of left-handed and right-handed leptons and quarks, the scalar slep- 
tons and squarks, and the corresponding antiparticles. Additionally, there are 
two Higgs doublets, the fermionic higgsinos and their antiparticles (see Ta- 
ble 2.2). Here, i = 1, 2, 3 is the generation index of (s)leptons and (s)quarks 
and h = r, g, b is the color index of (s) quarks. The table gives the notation 
of the particles in Lagrangians and Feynman rules, and their transformation 
properties under the Standard Model gauge groups. 

The chiral matter fermions carry spin 1/2 and their bosonic partners 
are scalars, while the Higgs bosons are complex scalars and the higgsinos 
are spin-1/2 fermions. The enlarged Higgs sector is needed to guarantee the 
cancellation of anomalies from the introduction of the higgsino superpartners. 
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Name Scalars 0* Fermions Xl {SU{3)c , SU{2)l) 

Sleptons, leptons D = (j}^,^ U = (^^^,^ (1, 2)_i 

E- = g^*^ E- = e^^^ (1,1)^1 



Squarks, quarks Ql = (^|' J Ql = i^J^'^J (3, 2)^. 

t^r = «Kh (3, 

^r = 4:/. ^^ = ^5. (3, 

Higgs, higgsinos = (^'] Hd = (1, 2)_i 



H+\ ~ fH. 



2 



Table 2.2: Chiral supermultiplets of the MSSM. 

In addition, the two Higgs doublets are needed to generate masses for 'up'- 
and 'down'-type quarks as well as charged leptons after electroweak symmetry 
breaking. 

Supersymmetric Lagrangians are determined by the superpotential W, 
which is a function of the supermultiplets.^ In the MSSM with conserved 
i?-parity the superpotential is given by [14] 

W4issM = fJ-HuHd + X^jHdLiEj + XfjHdQiDj — X^jHuQiUj . (2.2) 

The hats denote superfields and we sum over generation indices i, j = 1,2,3 
and the suppressed gauge indices, /i is the supersymmetric Higgs mass pa- 
rameter and A^-'^''^ denote the quark and charged lepton Yukawa coupling 
matrices. 



i2-Parity 

As a consequence of observed B — L conservation one usually introduces 
an additional parity called i?-parity in the MSSM. This parity is added by 
hand in order to forbid processes that lead to proton decay. Therefore, a new 
multiplicative quantum number 

= (_l)3(B-L)+25 (2.3) 

^We do not derive the MSSM Lagrangian here. Thus, we do not introduce the superfield 
formahsm and two component notation even if we write some expressions in this form. 
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is assigned to all the MSSM particles with spin S. Standard Model parti- 
cles have even {Rp = +1) and their supersymmetric partners odd -R-parity 
{Rp = —1). This imposes that supersymmetric particles can only be pro- 
duced pairwise at colliders and that they cannot decay into Standard Model 
particles only. 

Thus, the lightest supersymmetric particle would be absolutely stable. 
Moreover, to be consistent with cosmology, the LSP must be electromagnet i- 
cally and color neutral and therefore provides a natural candidate for particle 
dark matter. 

On the other hand, the most general MSSM superpotential contains also 
-R-parity violating terms [14]: 

Wf,^ = fi^HM + \\.,kUL,El + K.kUQM + ^K^.U^bpi. (2.4) 

As in the case of conserved i?-parity, summation over the generation indices 
z, j, = 1, 2, 3 and the suppressed gauge indices is assumed. The i?-parity 
breaking mass parameters /ij, and the trilinear Yukawa couplings Xijk and 
Kjk violate lepton number, while the couplings Af^-^. violate baryon number. 
There is no symmetry of the theory that forbids these terms. Experimentally, 
however, the couplings of these interactions are very restricted: The proton 
lifetime of Tp > 2.1 x 10^^ years [17] demands, for instance, that either the 
L violating or the B violating couplings vanish, or that all the couplings are 
extremely suppressed. In addition, the requirement that an existing baryon 
asymmetry in the early universe is not erased before the electroweak phase 
transition implies that Ajj^, A'^^-^ < 10"^ [33, 34, 35]. 

The distinction between Higgs and matter supermultiplets is lost in im- 
parity violating models. Therefore, i?-parity violation permits the mixing of 
sleptons and Higgs bosons, the mixing of neutrinos and neutralinos, and the 
mixing of charged leptons and charginos. In this case, the LSP is not stable 
anymore. This is a crucial fact for the decays analyzed in this work. 

Supersymmetry Breaking 

As mentioned before, supersymmetry must be spontaneously broken to ex- 
plain the mass differences of Standard Model particles and their superpart- 
ners. In order to solve the hierarchy problem, the SUSY breaking must be 
soft, i.e. the SUSY breaking parameters must not reintroduce quadratic di- 
vergences. Moreover, these parameters have to be at the TeV scale. 

There are no realistic models of spontaneously broken supersymmetry 
where the SUSY breaking arises from the particle interactions of the MSSM. 
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Usually one assumes a hidden sector, consisting of particles that are neu- 
tral with respect to the Standard Model gauge groups, and a visible sector 
that contains the MSSM particles. SUSY breaking is assumed to occur in 
the hidden sector and to be mediated to the MSSM by some mechanism. 
Supersymmetry is broken when one of the hidden sector fields obtains a non- 
vanishing F-term vacuum expectation value (VEV), labeled as {F). 

Since the SUSY generators are fermionic, the breaking of supersymmetry 
generates a massless goldstone fermion, the goldstino, analogous to the mass- 
less goldstonc boson of electroweak symmetry breaking. This is a problem of 
global supersymmetry since no massless fermion has been observed. 

The arising soft SUSY breaking parameters are, for instance, the mass 
parameters Mi, M2 and M3 of the electroweak and strong gauginos, respec- 
tively. Further parameters are, among others, the soft masses of the squarks 
and sleptons, and the soft masses rnfj^ and of the Higgs doublets. With 
broken i?-parity numerous additional soft terms arise, for instance the Higgs- 
slepton mixing parameters and m\.fj^. 

Electroweak Symmetry Breaking 

As we have stated before, in the MSSM there are two complex Higgs doublets. 
Electroweak symmetry is broken down to electromagnetism, 

SU{2)L^U{l)Y^U{l)em. 

dynamically through radiative corrections to the soft Higgs masses and 
mu^ - The neutral Higgs fields then acquire vacuum expectation values {H^) = 
Vu and (H^) = Va. The ratio of the Higgs VEVs is usually denoted as 

tan^ = — . (2.5) 

Vd 

The VEVs of the Higgs doublets are related to the Standard Model Higgs 
VEV V ~17A GeV in the following way: 

2 2,2 

Va = vsm/3, (2.6) 
Vd = V cos /3 . 

According to electroweak symmetry breaking through the Higgs mechanism, 
the electroweak gauge bosons absorb the three massless degrees of freedom 
of the two complex Higgs doublets: The Goldstone bosons G° and G^. These 
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become the longitudinal modes of the massive and vector bosons. 
The new mass eigenstates are given by 



Zl) \-^xxyew cosOwJ 
1 



(2.7) 

= ^ {W^ T iWl) , (2.8) 



where the weak mixing angle 9\y is defined as 



sm tJw 



cos 9 



w 



9 



(2.9) 



with the SU{2)l gauge coupling g and the U{1)y gauge coupling g'. The 
photon is massless, while the other electroweak gauge bosons acquire masses 

gv 

mw = -7= , 

^ gv (2-10) 

mz = — p . 

v2cos 6w 

The Higgs gauge eigenstates can be expressed in terms of the Higgs mass 
eigenstates: The neutral gauge eigenstates are decomposed as 

Hy - [vj + V2 [s. cj [hJ + V2 sj ^ ^ ^ 

and the charged gauge eigenstates read 



(2.12) 



Here and in the following parts we use the abbreviations Sq, = sin a and 
Ca = cosa for the mixing angles. In the tree-level approximation one has 
Po — /3± — f3 and the masses of the Higgs mass eigenstates are given by 

^h,H = \ {^\^ + ^1 T \j (m^o - + 4 m|m^o sin^ 2^^ , 

2 o I |2 , 2 , 2 (2-13) 

2 2,2 

m^± = m^o + 'm^^ . 
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In this case, the mixing angle a is determined by the conditions 

sin 2a ^ _ / m\^ml \ tan 2a ^ / m^p +m| \ 

sin 2/3 \^\~'^\) tan 2/3 \m^o — ""^1/ ' 

and is usually chosen to be negative. 

In the MSSM decoupling limit, i.e. for a large /i parameter, the particles 
i7, and are very heavy and decouple from low-energy experiments. In 
this case, using the relations (2.14), the mixing angle becomes a ~ /3 — 7r/2 
and the lightest Higgs boson h obtains the couplings of the ordinary Standard 
Model Higgs boson. We will use this framework later in this analysis. 



Neutralinos and Charginos 



Gauginos and higgsinos mix with each other due to electroweak symmetry 
breaking. The neutral gauginos and the neutral higgsinos combine to form the 
four mass eigenstates called neutralinos Xa , while the charged gauginos and 
the charged higgsinos mix to form the two mass eigenstates called charginos 
Xa- As mentioned before, the lightest neutralino Xi is a thoroughly studied 
candidate for the cold dark matter in MSSM models. 

In the gauge eigenstate basis '?/'° = (-B, W'^ 
part of the Lagrangian is given by 



ii^^ the neutralino mass 



1 



=5^cutralino mass 

where the neutralino mass matrix is symmetric and reads 



(2.15) 



( Ml 




TV 



9 Vd 

g 

V2 





M2 

gjJd 
V2 



9 Vd 


9 Vu 






9'"d 



















(2.16) 



Using the relations (2.6), (2.9) and (2.10), the mass matrix can be rewritten 
in the following way: 



N 



( Ml 


\ nizSwSfi 





M2 

mzCwC/B 
-mzCwS/B 



-mzswcp 
nizcwcis 


-/i 



mzswsi3 \ 
-mzcwsp 

/ 



(2.17) 



where we used the abbreviations sw = sin6'iy and cw = cos9w for the 
weak mixing angle. Using an orthogonal matrix that includes equation (2.7), 
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we can change the basis from the gauge eigenstates to the supersymmetric 
partners of the massive gauge bosons ip'^ = (7, 



Cw 



0\ 



Sw 

Sw Cw ^ ^ 

10 

\ 1/ 



R 



(2.18) 



The neutrahno mass term can then be rewritten in the form 

1 
2 

with 



-^^"^^M'^^"^ + h.c. 



(2.19) 



M'^ = RMnR^ 



2 



V 



2 

w 



(M2 - Ml) swcw 



(M2 - Ml) swCw Mislr + Msc^ mzcp 

mzC/B 
-mzs/s -fji 



\ 

-mzSi3 

/ 



(2.20) 



Similarly, the basis can be changed to the mass eigenstates, which are the neu- 
tralinos. Symmetric matrices can always be diagonalized by a unitary matrix 
S and its transposed matrix. Since we are interested in positive eigenvalues, 
we have to add an additional diagonal unitary phase matrix P (for the prob- 
lem of negative eigenvalues see for instance [36]). Using these matrices, the 
transformation into the neutralino basis reads 



\xV 



PS 



/7\ 



(2.21) 



so that the neutralino mass matrix is diagonalized in the following way: 



P*S*M'j^S^P^ 











\ 

























I 








(2.22) 



Similarly, for the chargino mass term in the Lagrangian we have in the 
gauge eigenstate basis — '"^^^^ ^i/- tj-\t 



1 



chargino mass 



(2.23) 
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where the chargino mass matrix Mc can be written in 2 x 2 block form 





^-=\X J ' (2.24) 
with 

For the last expression we used the relations (2.6) and (2.10). We can change 
the basis to the mass eigenstates using two unitary matrices U and V: 

(I) -(I)- 

The matrix X can be diagonalized by U and V in the following way: 

U*XV^ = VX^U^ =\^^ ° I, (2.27) 
where the masses are the positive roots of the eigenvalues of X^X, since 

VX^XV^ = U*XX^U^ = \ 2 , (2.28) 

and are given by 

ml^^ = U\M2\^ + \l^f + 2mlr 



T V(l^2| +2m2^,) -4|/xM2-m2^sin2/3|' 



(2.29) 



In models where the gaugino masses unify at the gauge coupling unifica- 
tion scale Mjj, one has the tree- level relation [25] 

5 

Ml ^ - tan^ ewM2 ^ 0.5 M2 (2.30) 

3 

between the bino and the wino mass parameter at the electroweak scale. In 
this work we will use the GUT relation 

M2~1.9Mi, (2.31) 

obtained from the numerically computed renormalization group evolution of 
the gaugino masses down to the electroweak scale [5]. 
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Name 



Bosons Fermions {SU{3)c , SU{2)l) 



Y 



Graviton, gravitino 



9^ll^ i^i^ (1, l)o 



Table 2.3: Gravity supermultiplet. 



2.2 Supergravity 

If supersymmetry is promoted to a local symmetry, i.e. the parameter in 
SUSY transformations becomes coordinate-dependent, the theory must in- 
corporate gravity. This is because in order to achieve invariance under local 
SUSY transformations, one has to add a new supermultiplet to the theory: 
The gravity supermultiplet, which consists of the spin-2 graviton and the 
spin-3/2 gravitino (see Table 2.3). The resulting locally supersymmetric the- 
ory is called supergravity (SUGRA). The gravitino, as well as the graviton, 
is neutral with respect to the SM gauge groups and has vanishing mass (and 
therefore only two transverse helicity states) in the case where supergravity 
is unbroken. 

Since in supergravity there appear couplings with negative mass dimen- 
sion, the theory is nonrenormalizable. On the other hand, we assume that 
SUGRA is an appropriate low-energy approximation of a more general the- 
ory. In the flat limit, i.e. Mpi — )■ oo, the renormalizability is restored. In this 
work, however, we are interested in couplings of the gravitino to MSSM par- 
ticles that are suppressed by the Planck mass. Therefore, we do not work in 
the flat limit, but take into account only tree- level interactions. 

Supergravity Breaking and the Super-Higgs Mechanism 

Analogous to the Higgs mechanism of electroweak symmetry breaking, in 
supergravity there exists a super-Higgs mechanism of supergravity breaking. 
The massless Goldstone fermion of supersymmetry breaking, the goldstino, 
is absorbed by the gravitino, which thereby acquires its longitudinal (helicity 
±1/2) components and becomes massive. 

This is because in spontaneously broken supergravity the Lagrangian con- 
tains a gravitino-goldstino mixing mass term. Invariance of the gravitino and 
the goldstino fields under local supersymmetry transformations demands a 
redefinition of the fields. The redefined gravitino is a linear combination of 
the gravitino and the goldstino, and therefore gets all four helicity states. In 
this case the gravitino mass becomes 



(2.32) 



m3/2 



V3Mp ' 
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where {F) is the F-term VEV responsible for the spontaneous breaking of 
supergravity. 

The value of the gravitino mass depends on the particular SUSY breaking 
scheme and can range from the eV scale to scales beyond TeV. For instance, 
in gauge-mediated SUSY breaking the gravitino mass is usually much less 
than 1 GeV, while for gravity mediation a mass in the GeV to TeV region is 
expected [37]. Regardless of the specific SUSY breaking mechanism, we take 
the gravitino mass as a variable parameter that is expected to be about (9(10- 
100) GeV in order to account for the dark matter density of the universe (cf. 
Section 2.3). 



2.3 The Gravitino 

In the following we will discuss the classical gravitino field and its conse- 
quences for the standard cosmological scenario. 

2.3.1 The Free Massive Gravitino Field 

The massive gravitino is described by the following Lagrangian [38]: 

^ = '1 ^^""""i^pl^'ludp'^. - \ my2^^ [7^ Y] + ^^nt . (2.33) 

For the free gravitino field we discard the interaction part of the Lagrangian. 
Application of the Euler-Lagrange equations to this Lagrangian results in the 
equations of motion for the free massive gravitino field: 

= — Oy 



di)^ didyijj^) ^2.34) 
= -I e^'""'l'ludpi^. - \ ^3/2 [7^ 7l = . 

From this expression a simpler set of equations of motion can be derived. To 
achieve this, we apply either or the identity 

^^e'-'^p^ = ^^.y^. _ ^p.^u ^ gu.y _ gup^.^ ^ (2.35) 
that is derived in Appendix B, to equation (2.34). We obtain 

^3/2 (^7^^^^. - 7'^'^^.) =0, (2.36) 

I {^^$Y^u - rd'^y + $V - Yd^ijy) + m3/2 (7^7'V^. - V) = ^ ■ (2-37) 
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Application of 7^ to (2.37) then results in 

I [h'^u - I'Hu) + 3 m3/27>. = . (2.38) 

From equations (2.36) - (2.38) the Rarita-Schwinger equations [39] for the 
massive gravitino field can be derived: 

7'^^^(x) = 0, (2.39) 
{%$ - ms/s) ^^(x) = . (2.40) 

These equations yield the further constraint 

<9'^V^^(x) = 0. (2.41) 
From (2.39) - (2.41) we can also derive the adjoint equations: 

V';.(x)7^ = 0, 

idvi^Ml" + ^^I2^M = , (2.42) 
9^V^^(x) = 0. 

In order to find the solutions to the Rarita-Schwinger equations, we go to 
momentum space and note that similar to the case of Dirac spinors there are 
positive and negative frequency solutions:^ 

M^) = V^r(P) e"'"" ^,{x) = r,'{v) e^^■^ s = ±1 ±^ , (2.43) 

where the mode functions ip'^ and ip~ have to obey 

7'^V^+^(p) = 0, 7/^V^;^(p) = 0, (2.44) 

{p-m3/2)^t'{p) = 0, and {p + my2)^;'ip)=0, (2.45) 
P^i^t'{p) = p/^^;^(p)=0. (2.46) 

Following [40] the mode functions can be constructed using the tensor product 
of the familiar Dirac spinors u and v of spin-1/2 particles, and the polarization 
vector of a massive spin-1 particle. The result is 



l,s])u"^{p)e'^{p) 
l,s])v-{p)e'^{p) 



(2.47) 



^ As stated in [30] , all gravitino processes can be described using only one mode function, 
since the gravitino is a Majorana particle. However, in order to maintain the freedom to 
choose the direction of the continuous fermion flow in the Feynman rules (see Appendix C), 
we have to discuss both mode functions. 
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1712 = +1 ^2 = 1712 = —1 



mi = +i 1 

mi = -| y/ys v/273 1 

Table 2.4: Clebsch-Gordan coefficients for ji = ^, ^'2 = 1 and J = |. Figures 
taken from [17]. 

where ((ji, mi) (j2, ""^2)1 {J, M)) are Clebsch-Gordan coefficients with M = 
mi + m2. The values of the coefficients in equation (2.47) are given in Ta- 
ble 2.4. 

Using this form of the mode functions, and the normalization of the Dirac 
spinors (B.25) and the polarization vector (B.31), we can derive the normal- 
ization of the gravitino mode functions. The result is 

^P^^ip)r''''ip) = -2m,/2 5'', 
^-^{p)^-^'^{p)=2m,/25'''. 

In the calculation of unpolarized decay rates we encounter a summation over 
the helicity states of the gravitino field. The gravitino polarization tensor is 
defined as 

s 

where the sum is performed over the four gravitino polarizations s = ±|, ±i. 
Using the normalization of the mode functions (2.48) we obtain 

p^xip)p^Mp) = J2'^i'ip)^rip)i^^''\p)^^''ip) 

(2.50) 

= T2m3/2Pi(p). 

Thus the polarization tensor is a projector. We find that the polarization 
tensors for a gravitino with four-momentum p are given by 

P^P) = -if + ^3/2) |n^.(p) - I U^Ap)U,,{p)ri'^ (2.51) 
for the positive frequency mode functions and 

P;:uip) = '{f- ^3/2) - \ n^.(p)n.A(p)7V} (2.52) 

for the negative frequency mode functions. In the above expressions we use 

n,.(p)= U,.-^). (2.53) 



"^3/2^ 
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For a derivation of the polarization tensors see Appendix B. 

Since the gravitino field is a solution of the Rarita-Schwinger equations 
of motion, the polarization tensor obeys the following constraints 

7^Pi(p) = 0, P^^(p)Y = 0, 

p^P%{p) = ^, and P±(p)p- = 0, (2.54) 

T m3/2) PtM = PtM {f T m,„) = . 



2.3.2 Gravitino Interactions 

In the previous section we discussed the free gravitino field. Now we want to 
add interactions between the gravitino and the MSSM fields. The gravitino 
interaction part of the Lagrangian reads [29, 38] 

i 



V2Mpi ^2.55) 



-2\ 



We immediately see that all the gravitino interactions are suppressed by the 
Planck scale. The covariant derivative of scalar fields is given in the MSSM 
as [29] 

3 

D,<P. = d,<P, + I J2 ^?«4"^ '^^i'^i'^^- (2-56) 

a=l 

and the field strength tensor for the gauge bosons reads 

« = d^Al"'^ " - ) " - '^^'^^(f ) ''Al"'^ ^ (2.57) 
T^^ij, a = 1, 2, 3 are the generators of the Standard Model gauge groups 



(2.58) 



where Yi is the hypercharge. The Pauli sigma matrices cTq are given in (B.7) 
and Xa are the eight Gell-Mann matrices which will not be needed in this 
work. are the totally antisymmetric structure constants of the corre- 

sponding gauge group. 



7^(1) 

a,ij 




a, ij 


1 


rp{3) 

a,ij 





32 



2.3. The Gravitino 



The Feynman rules can be extracted from the interaction Lagrangian 
in the usual way. Since the gravitino and the gauginos are Majorana fields 
for which exist Wick contractions different from those of Dirac fermions, 
amplitudes will contain charge conjugation matrices and there may arise 
ambiguities concerning the relative sign of interfering diagrams. 

Therefore we use a method that introduces a continuous fermion flow [41]. 
The direction of this fermion flow in a process can be chosen arbitrarily, if the 
corresponding Feynman rules are used. Amplitudes are then written down in 
the direction opposite to the continuous fermion flow. This method avoids 
the appearance of charge conjugation matrices and gives the correct relative 
signs of different diagrams contributing to a single process. 

The complete set of Feynman rules that are relevant for this work is 
provided in Appendix C. 



2.3.3 Gravitino Cosmology 

The addition of the gravitino to the particle spectrum can lead to several 
cosmological problems. We want to discuss shortly two of the main problems 
and see what restrictions for gravitino theories arise therefrom. 



Thermal Gravitino Production 

If gravitinos are in thermal equilibrium in the early universe, their relic den- 
sity leads to overclosure of the universe. That means they contribute a density 
^3/2 = f?3/2/f?c > 1- This problem can be circumvented in an inflationary uni- 
verse, since during inflation any initial abundance of gravitinos is diluted due 
to the exponential expansion of the universe. Gravitinos can then be repro- 
duced by scattering processes in the thermal plasma after the universe has 
been reheated (cf. Section 1.1). Thermal gravitino production is enhanced 
due to the contribution from the less suppressed interactions of the goldstino 
component. The gravitino relic density in that case becomes [38] 

^ ,9 f Tr \ /100GeV\ / ms \2 , , 

where Tr is the reheating temperature of the universe after inflation and rrig 
is the gluino mass. For reasonable values of the gluino mass and a gravitino 
mass in the (9(100) GeV range, thermally produced gravitinos can amount 
to the observed dark matter density, i.e. ^23/2 — ^dm- This constrains the 
reheating temperature to be Tr ^ (9(10^°) GeV, which is compatible with 
the constraint Tr > 10^ GeV from thermal leptogenesis (cf. Section 1.1). 
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Impact on Big Bang Nucleosynthesis 

Long-lived supersymmetric particles can affect the successful predictions of 
BBN if their lifetimes are longer than ~ 1 s, i.e. if they decay during or after 
primordial nucleosynthesis. The particles produced in these late decays may 
induce hadronic and electromagnetic showers and thereby lead to hadro- and 
photo-dissociation of light elements. 

The gravitino — in cases where it is not the LSP — decays into lighter par- 
ticles with a lifetime [42] 



The late decay of the gravitino then implies an upper bound on the reheat- 
ing temperature to ensure that the BBN predictions are not significantly 
altered [43]: Tr < 10^ GeV. This low value for the reheating temperature is, 
however, not compatible with the value favored by thermal leptogenesis and 
the paradigm of gravitino dark matter. 

On the other hand, if the gravitino is the LSP and stable, the lifetime 
of the next-to-lightest supersymmetric particle decaying into a gravitino and 
Standard Model particles is given by [25] 



the predictions of standard BBN. For instance, the hadronic decays of a neu- 
tralino NLSP can dissociate the primordial generated light elements, whereas 
a long-lived stau NLSP can form a bound state with ^He and catalyze the 
production of ^Li. This can lead to an overproduction of ^Li by a factor 300- 
600 [44]. Other possible NLSPs like sneutrinos or stops do not substantially 
affect the standard BBN predictions [45, 46]. 

However, there is a more general solution to this problem. The introduc- 
tion of a small i?-parity violation causes the NLSP to decay into Standard 
Model particles before the onset of BBN. Due to the double suppression of 
the gravitino couplings to SM particles by the Planck mass and the small 
i?-parity violation, the gravitino remains very long-lived and therefore a vi- 
able candidate for the cold dark matter [47]. This paradigm we will further 
discuss in the following chapter. 




(2.60) 




(2.61) 



Thus, the NLSP is present during or after BBN and its late decays may spoil 
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Gravitino Decays via Bilinear 
ii-Parity Breaking 

In this chapter we will shortly introduce bilinear i?-parity breaking models 
of the type discussed thoroughly in [4] and then discuss in detail the decay 
channels of the LSP gravitino via such i?-parity violating couplings. 

3.1 Bilinear iZ- Parity Breaking 

In models with bilinear i?-parity breaking, the MSSM superpotential contains 
the additional term 

W^^ = li,UHu, (3.1) 

where /Xj are i?-parity violating higgsino-lepton mixing masses and summa- 
tion over the generation indices i = 1, 2, 3 is assumed (cf. Section 2.1). 

From equation (2.3) we see that i?-parity breaking implies B — L breaking. 
Thus, it can be shown that /Xj is suppressed compared to the supersymmetric 
higgsino mass parameter by ~ f^_^/Mp[, with vb-l being the scale 
oi B — L breaking. The i?-parity violating bilinear coupling Hi can be rotated 
away by a redefinition of the Higgs and lepton superfields: 

Ha = H',-^L\ and U = L[ + ^H',. (3.2) 

This mixing of the Higgs and lepton superfields induces trilinear i?-parity 
breaking Yukawa couplings 

Kjk = —Kj and A- -fc = — Af- , (3.3) 
/i fj, 

that are suppressed by 0{fii/ fi). The baryon number violating couplings A"^-;. 
are suppressed at higher order in this model. Thus, the model is compatible 
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with the bound from the proton hfetime. In order to be also compatible 
with cosmology, there are additional constraints on the Yukawa couplings: As 
already mentioned in Section 2.1, the viability of baryogenesis via leptogenesis 
requires A, A' < 10~^. Otherwise the generated baryon asymmetry is washed 
out before the electroweak phase transition. On the other hand, the couplings 
must be large enough to ensure that the NLSP decays before the onset of 
BEN. A sufficiently short lifetime can be achieved with A, A' > IQ-^^ [4]. 

Apart from the supersymmetric term discussed above, the corresponding 
soft bilinear supersymmetry and -R-parity breaking terms Bi and Tn\_^^ arise 
in the Lagrangian: 

-^sojt = m%jHd\'^ + m]jJ\Hu\^ + m\\^Li\^ 

+ (^BHdHu + BiLiHu + ml^H^LiH*^ + h.c^ + . . . . ^^'^^ 

Minimization of the scalar potential yields besides the Higgs VEVs a non- 
vanishing vacuum expectation value along the sneutrino field direction: 

{h) - • (3.5) 

This sneutrino VEV explicitly breaks lepton number and generates not only 
one neutrino mass, but also nonvanishing mixing between neutralinos and 
neutrinos, as well as between charginos and charged leptons [14]. In addition, 
the sneutrino VEV can be used to parameterize the Higgs-slepton mixing. 
For later convenience we thus rewrite the relation for the sneutrino VEV in 
the following form: 



5isin/3 + mi^j;^^cos/3 ~m?^-^. (3.6) 



The above mentioned mixings are responsible for the two-body decays 
of the gravitino into gauge/Higgs boson and neutrino, which are the main 
source of neutrino fiux in our scenario. These decays are also possible at the 
one-loop level if only trilinear i?-parity breaking terms are considered [48], 
but we do not consider that case here. Since the neutralino-neutrino mixing 
takes place along the zino component, the branching ratios into the different 
gauge/Higgs boson channels are fixed by the neutralino mixing matrix once 
the gravitino mass is specified. 

In this analysis we assume that the sneutrino acquires a VEV only along 
the direction, since we expect the -R-parity violating couplings to be largest 
for the third generation. The produced neutrinos then have tau fiavor. How- 
ever, we will see later that this is not a crucial assumption, since neutrino 
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oscillations change any pure neutrino flavor into a mixed state. In particu- 
lar, due to maximal atmospheric mixing, the flux of tau and muon neutrinos 
turns out to be identical (see Chapter 4). 

3.2 Gravitino Decay Channels 

In this section we compute the relevant decay channels of the LSP gravitino 
in the above introduced model, where the neutralino-neutrino mixing, as well 
as the chargino-charged lepton mixing, is generated through a nonvanishing 
tau sneutrino VEV {vt). The primary decay channels of the gravitino into 
neutrinos and other Standard Model particles are 

-03/2 -^IT^r, 

^3/2 ^ W^r^, 

The flrst decay mode is practically always allowed, since the tau neutrino 
mass is very small. On the other hand, the decay modes into W^t^ and 
Z^Ut- are only accessible for a gravitino heavier than the weak gauge bosons, 
while the channel into h requires a gravitino mass above the threshold for 
the production of the lightest Higgs boson. 

The total decay width of the gravitino determines the gravitino lifetime 
according to 

n/2 = ^ ■ (3.8) 

J- 3/2 

However, the partial decay widths of the gravitino will contain the sneutrino 
VEV from bilinear R-parity breaking whose value is not known. Thus, we 
will not be able to flx the absolute lifetime of the gravitino. Nevertheless, 
from the relative decay widths of the different decay channels we are able to 
determine the branching ratios for these channels. This is a crucial input for 
the prediction of spectra of gravitino decay products. 

3.2.1 Calculation of Decay Widths 

The gravitino decay widths can be computed from the Feynman rules given 
in Appendix C. Here we will only present the detailed calculation for the 
decay width of the process 1P3/2 1 ■ The calculation of the other decay 
widths is given in Appendix C. 
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V^3/2 -^IJ^T 

At tree level only one Feynman diagram contributes to the decay of the LSP 
gravitino into the photon and the tau neutrino: 




There is no direct coupling of the photino to the tau neutrino, since it is elec- 
tromagnetically neutral. However, the photino mixes with the zino through 
the neutralino mass eigenstates and the zino couples to the tau neutrino via 
the tau sneutrino VEV. Using the rotation of the neutral gauge eigenstates 
into the neutralinos (cf. Section 2.1) 

4 4 

7 = E and Z'=Y. ^h^Lxl , (3.9) 

i,a=l i,a=l 

and the relation between the gauge eigenstates and mass eigenstates of the 
electroweak gauge bosons (2.7) together with the generators of the corre- 
sponding gauge groups (2.58), we can write the amplitude of this process 
as 

iM =- u''{q) (z>^) (fi'-y^ cos dw - sin (^w^ Pr 




PI 



To approximate the gaugino propagator we used the fact that = ml <^ 
m|o • Although the tau neutrino is no mass eigenstate we regard it as the final 
state here. Its mass is depicted as m,^ and is negligible in this calculation. In 
addition, we introduced gz = g/ cos9\y. In the next step we will introduce 
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= iyr) /v and the photino-zino mixing parameter 



Since we are looking at the decay of a gravitino and therefore work in the 
gravitino rest frame, the gravitino has no defined hehcity. Therefore, we have 
to average over the initial gravitino spin states. And since we do not want 
to measure the polarization of the decay products we sum over their helicity 
states. Then we can use the polarization sums for gravitinos (2.51), fermions 
(B.26) and polarization vectors (B.32), and the squared amplitude becomes 



256 I 



■ Tr [(^ + m,) (1 + 7^) ^ Y] P^P) [7^ ^] 7^ (l " 7^)] 

" tS^, I^^^I' + ^'^'^^^^ f^^' ^] 



In the second step we used the projector property of the chirality operators 
(B.16) and in the last step we already replaced squared four-momenta by the 
corresponding squared particle masses. The trace could be further simplified 
and calculated by hand exploiting the identities (B.40) and (B.41) for the 
gamma matrix structure and the constraints (2.54) for the polarization ten- 
sor. However, here we used for the calculation of the trace the Mathematica 
package Feyncalc [49]. 

According to equation (C.3), the further scalar products of the four- 
momenta for this process are given by 

{p-k)= "\ '^ =(fc-g), 
m2 + ml 

{p-q) = ^ • 

Thus, the squared amplitude becomes 

2 H/2 - (3 "^3/2 + 2 
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Using equations (C.7) and (C.8) we finally obtain the decay width: 

3 



r (^^3/2 7 I^t) 



^3/2 |2 I _ 



1 ml 
1 + ^ 



'3/2 



3 mgyg 



(3.11) 



or if we neglect the contributions of the small neutrino mass 



(3,12) 



The conjugate process ip3/2 l^r, that produces tau antineutrinos, has the 
same decay width: 




The decay widths of the other channels are given by^ 



mz 



(3.13) 



3 ms/2 6 



r {i'3/2 h Ur) 



"^3/2^4 
2 Ph 
PI 



3847r M, 



nit 



TTit, — mr 



+ sin (3 Ugo z + cos f3 Ufro v 

d 



(3.14) 



where the zino-zino, wino-wino and higgsino-zino mixing parameters are 



^ These results do not exactly coincide with those presented in [7]: For the interference 
terms proportional to jx we find a negative sign and a larger coefficient. 
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given as 



Uzz = ^zY^ 



4 Q* c 



a=l 



2 ^ m-± 



Q. = l Aq, 



and the kinematic functions fx-, jx and /ix are given by 



/3x = l- 
/x = l + 
Jx = 1 + 



m 



X 



2 

"^3/2 



2 1 



3 m3^2 3 777-3^2 

1 



(3.16) 



2 7773^2 
,2 



/l V = 1 + 10 ^ + 



7713^ rrix 



2 4 ■ 

"^3/2 "^3/2 

The decay widths for the conjugate processes ip3/2 Z^i^r, V'3/2 W^r^ 
and 4'3/2 h u^. are equaL 

The Higgs decay channel and the contribution from the non-abehan 4- 
vertex to the partial widths of the gravitino decay into the weak gauge bosons 
was previously neglected in [5]. Therefore the contribution of the photon 
channel has been overestimated. On the other hand, it seems that the photon 
channel has been underestimated in [7] due to the different result for the 
interference term in the and W decay channels. 

Even if we cannot determine the absolute value for the gravitino lifetime 
due to the unknown value of the tau sneutrino VEV, we are still able to 
determine the branching ratios of the different decay channels. They are 
defined as the ratio of the partial decay width and the total decay width: 

rtot = Vr(V^3/2->X) and BR(X) = ^^"^'(f "^^^ . (3.17) 
^ -L tot 

Using the same set of SUSY parameters as [7] — i.e. we work in the MSSM 
decoupling limit (see Section 2.1, we use = 10 TeV for the numerical calcu- 
lation) and use tan/3 = 10, nih = 115 GeV, mj> = 27773/2, Mi = 1.57773/2 and 
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Figure 3.1: Branching ratios of the different gravitino decay channels as a 
function of the gravitino mass in the case of the MSSM decouphng hmit. 

the GUT relation (2.31) between Mi and M2 at the electroweak scale — we 
can calculate the branching ratios for the gravitino decay channels. For the 
weak mixing angle we use the ratio of the electroweak gauge boson masses 



The result for the branching ratios is given as a function of the gravitino mass 
in Figure 3.1 and for several specific masses in Table 3.1.^ As expected, the 
branching ratio for the decay into photon and tau neutrino drops drastically 
for gravitino masses above the threshold for the production of weak gauge 
bosons. This is explained by the i?-parity violating mixing of the zino and 
wino to the tau neutrino and the tau lepton, respectively. The photino, on 
the other hand, has no direct coupling and the decay width is therefore 
suppressed with respect to the weak gauge bosons. 

For the above results we used a numerical calculation of the mixing pa- 
rameters [/-^, Uzzi U^ri^ and ?7„o 7- Following [5], we can also give an 

u, d 

analytical approximation for some of these mixing parameters. In the MSSM 
decoupling limit, the photino, zino and wino are approximately mass eigen- 
states since the off-diagonal terms in the neutralino and chargino mixing 
matrices are subdominant (see Section 2.1). In this case, the zino-zino and 
wino-wino mixing parameters are simply given by the propagators of those 

^Due to the different result for the interference terms in equation (3.14), the photon 
channel in our result is less suppressed than in [7]. 



cos6'- 
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"^3/2 








BR{h Pr) 


lOGeV 


100% 








85 GeV 


53% 


47% 






100 GeV 


7.5% 


83% 


9.0% 




150 GeV 


1.1% 


70% 


28% 


1.2% 


250 GeV 


0.34% 


60% 


29% 


11% 


1000 GeV 


0.021 % 


51% 


25% 


24% 



Table 3.1: Branching ratios of the different gravitino decay channels for sev- 
eral specific gravitino masses in the MSSM decoupling limit. 



particles, rescaled by the Z° and W boson masses, respectively: 
\Uzz\ - T7-> = n.. and \U^^\ ~ —- = — - . (3.19) 

As the photino has no direct coupling to the neutrino, the mixing parameter 
is approximated by the product of the photino propagator, the photino-zino 
mixing and the zino propagator, rescaled by the Z° boson mass: 

I I "^^Knz _ mz (M, - M,) swcw . . 

These approximations differ by less than 10 % from the numerical calculation. 
The ratios of the mixing parameters are then given approximately as 

|f^^zh|f^zzh|f^wl ^1:3.2:2.6, (3.21) 

where we only had to imply the GUT relation (2.31). For the higgsino-zino 
mixing parameters this method cannot be used since and are no mass 
eigenstates. However, in the MSSM decoupling limit the higgsino-zino mixing 
parameters are suppressed due to the large /i parameter. In that case, the 
diagram with Higgs-slepton mixing dominates the decay channel. For tau 
sneutrino masses that are not degenerate with the lightest Higgs mass we 
have m?^/(m?^ — m^) = 0(1), and since the tau sneutrino has to be heavier 
than the LSP gravitino, we conclude that there is no strong dependence of 
the Higgs channel on iriy^. 

We see that the branching ratios given in Figure 3.1 and Table 3.1 are 
not only valid for the specific choice of parameters listed above, but approxi- 
mately also for a large parameter space. In fact, the branching ratios for the 
electroweak gauge bosons dominantly depend only on the gravitino mass, 
while the Higgs channel additionally depends on the Higgs mass by means of 
the phase space factor. 
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Decay mode 


Branching Ratio 




e+ e 


(3.363 ± 0.003) % 






fo Odd _i_ n nn'7\ 0/ 
(^O.OOO ± U.UU/^ j 70 






(^o.o /^U ± U.UUo j /o 




invisible 


/'on nn _i_ n n^N 0/ 
(^zU.UU ± U.Uoj 70 




hadrons 


(by. 91 ± O.Odj 70 


w- 


e' Ve 


(10.75 ±0.13)% 






(10.57 ±0.15)% 






(11.25 ±0.20)% 




hadrons 


(67.60 ± 0.27) % 






(17.36 ± 0.05) % 






(17.85 ± 0.05) % 






(10.91 ±0.07)% 






(25.52 ±0.10)% 




TT^ Tl^ Tl^ Vj. 


(9.27 ±0.12)% 




TT^ 7r+ TT^ Z/,- 


(8.99 ± 0.06) % 



Table 3.2: Dominant decay modes of the Z and W bosons, and the r lep- 
ton. The conjugate decay modes have equal branching ratios. Figures taken 
from [17]. 

3.2.2 Fragmentation of the Z*^, W and h Bosons 

Except for the photons and neutrinos, the decay products of the above dis- 
cussed gravitino decay modes are not stable. Since we are interested in gravi- 
tino decays at astrophysical distances, we have to determine the spectra of 
stable particles at the end of the fragmentation processes of the W and 
h bosons and the decay of the r lepton. The stable final state particles are 
photons, electrons, protons, neutrinos and their corresponding antiparticles. 

The dominant decay modes and the corresponding branching ratios for 
the first fragmentation step are presented in Table 3.2 and in Figure 3.2. 
Since we work in the MSSM decoupling limit, the lightest Higgs boson h has 
the same interactions as the Standard Model Higgs boson. 

In fact, the complete calculation of fragmentation processes is very compli- 
cated. To obtain the decay products from the Z^v^, W^t^ and h channels 
and their energies, we therefore mimic the decay of a resting gravitino with 
the event generator PYTHIA 6.4 [51]. A description of the used PYTHIA 
program is presented in Appendix D. 

In order to infer the neutrino spectra from the neutrino energy distribu- 
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Decay 


mode 


Z/e + Ue 










16.2 


29.6 


1.2 






15.9 


28.9 


1.2 




h Ur 


29.5 


54.1 


1.2 



Table 3.3: Numbers of neutrinos and antineutrinos from the fragmentation 
of the different decay products simulated with PYTHIA. 



tion obtained from PYTHIA, we need a large number of generated events for 
the simulation. Thus we use 10^ events for most neutrino channels. For the 
strongly peaked tau neutrino channels of the gravitino decay into Z^Vj. and 
h Vr we use twice the number of events. 

The numbers of neutrinos of different flavors, generated in the fragmen- 
tation of the different decay products, are summarized in Table 3.3. We note 
that a large number of muon and electron neutrinos is produced in the frag- 
mentation. At lower energies their ratio is 2:1 as in the case of atmospheric 
neutrinos (cf. Chapter 4). 
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3.2.3 Neutrino Spectra from Gravitino Decays 

The injection spectrum of neutrinos from gravitino decays is composed of a 
series of contributions. The gravitino decay into a photon and a tau neutrino 
is a two-body decay with a very long hfetime and virtually vanishing daughter 
particle masses. Therefore, the resulting neutrino spectrum from this channel 
is just a line at half the gravitino mass: 



dE 



(V'3/2^7^r)^<^(£^-^). (3.22) 



By contrast, the decays into the heavy gauge bosons and the Higgs boson 
contribute a continuous spectrum. Anyway, the main features of the frag- 
mentation spectra are also lines. The decays of the gravitino into Z^Vt and 
h v-j. result in lines centered at 




E^^^'^ll--^]. (3.23) 



However, these lines are not monoenergetic, but have a shape described by a 
normalized Breit-Wigner profile: 

— ^ (V'3/2 ^ h = ^ 9 — — , 3.24 

where the widths are given in terms of the /h boson widths according to 



dE, 



dmzh 



rz,.= ^r^,.. (3.25) 



The leptonic decays W ^ Ivi, Z^ ^ viVi and h ^ Ivi also produce, in the 
rest frame of the decaying particle, monoenergetic neutrinos in all flavors. 
However, due to the boost of the bosons in different directions, the lines smear 
out almost completely in the Earth's rest frame, giving just an additional 
contribution to the continuous part of the spectrum. 

As already mentioned in the previous section, we extracted the continuous 
fragmentation spectra for the different neutrino flavors in the W^t^, Z^Vr 
and h channels from the PYTHIA simulation. We denote these spectra by 
dN^'^/dE, dN^^^/dE and dN^^^'/dE, respectively. Weighting the constituent 
channels with the corresponding branching ratios, the total spectra for the 
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Figure 3.3: Neutrino injection energy spectra from gravitino decay for the 



different flavors in the case of m 



3/2 



150 GeV. The dominant hne is 



clearly visible in the tau neutrino spectrum at an energy of 47 GeV, while the 
photon line is located at the end of the spectrum. Considering a Higgs mass 
of 115 GeV, the peak at ~ 30 GeV is strongly suppressed by the phase space 
factor. 



(3.26) 



different neutrino flavors are given by 

dN^ , , dN"^^ , n X dNf/ , , rfV^^ 

_p =BR(W--r-) + BR(ZV) + BR(A.,)-^ 

dN, , ^ . dN'^^ , , , dN^^ dNj"" 

-IW- - f^"'^^) ^ + BR (ZO.,) + BR („ 

+ BR(ZV)^ + BR(A.,)-^. 

The resulting energy spectra for the different neutrino flavors are shown in 
Figure 3.3. In this case, and in most parts of this analysis, we consider a gravi- 
tino mass of m3/2 = 150 GeV, motivated by the interpretation of anomalies 
in other cosmic ray channels as mentioned in the preface. Furthermore, at 
lower gravitino masses, the detection of any signal is much more difficult due 
to the lower neutrino yield and higher background fluxes. 

For this value of the gravitino mass, both the Z° and Higgs lines are visible 
in the spectrum in addition to the photon line. We see that the tau neutrino 
spectrum shows a very characteristic signature of two or three distinctive 
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peaks (three in the region above the hghtest Higgs threshold, as in this case) 
in addition to a continuum which is suppressed at low energies. From Figure 
3.1 we see that the dominant line from the decay into Z^Vj. always has a 
branching ratio larger than 25 %. The spectra of the other two flavors are very 
similar to each other, consisting only of a continuum contribution practically 
following a power law behavior oc E~'^ at energies below the sharp threshold 
at half the gravitino mass. 
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Chapter 4 



Neutrinos as a Probe of 
Gravitino Dark Matter 



In the previous chapter we have seen that the decay of the LSP gravitino 
in i?-parity breaking models produces stable Standard Model particles with 
characteristic energy spectra. The consequences of astrophysical gravitino 
dark matter decays have already been studied in several works: For instance, 
the gamma-ray signal from gravitino decays is studied in [5, 24], while the 
antimatter signals (i.e. fluxes of positrons and antiprotons) are discussed 
in [6, 7, 52]. 

In this work we consider the neutrino signal at the Earth arising from 
gravitino dark matter decays and discuss the implications for the detection 
of these neutrino fluxes. In addition, we will use the neutrino signal to place 
constraints on the gravitino parameter space. 

4.1 Neutrino Signals 
4.1.1 Neutrino Fluxes 

In this section we determine the neutrino flux at the Earth originating from 
gravitino decays. In this context the flux is defined as the number of particles 
per unit area and unit time. In fact, however, in most parts we will work with 
the differential flux per unit solid angle: 



There are two sources for a diffuse neutrino background: Gravitinos decaying 
in the dark matter halo of the Milky Way and those decaying in the halos of 
other galaxies at cosmological distances. 



_ d(f) 

^ dn 



dAdt dn 



dN 



(4.1) 
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Extragalactic Contribution Gravitino decays at cosmological distances 
give rise to a redshifted spectrum. A neutrino, generated in a gravitino decay 
at the comoving distance x{^) with energy Eq, is observed on Earth with 
an energy E = Eq/{1 + z) due to the redshift of relativistic particles in 
an expanding universe. The number of generated neutrinos per unit energy 
and time in a comoving volume element located at the comoving coordinates 
X = ix, 0, 0) is given by 

dEodt T3/2 dEo T3/2 dEo 

Here, the number density of gravitinos n3/2{x) = ^3/2Qc/^3/2 is assumed to 
be constant on cosmological scales and dN^/dEo is the neutrino energy spec- 
trum from a single gravitino decay (cf. previous chapter). Since the gravitino 
lifetime is much larger than the age of the universe, the decrease of 71^/2 
through decays can be neglected. 

An observer at the origin of the coordinate system receives a differential 
flux per unit energy of 

00 00 
dJeg _ f dx dNo{x) _ fl3/2Qc f dNi 



dE 



r dx dNoix) _ ^3/2Qc f dN, 

J A-Kx^dEdtdn 47rr3/2m3/2y d{E{l + z)) ^' ^ ' ' 




In order to perform the integration over the comoving distance, we have to use 
the dynamics of the universe presented in Section 1.1. Neutrinos propagate 
very similar to light since they are neutral and highly relativistic. Thus, 
we have for the line element according to the Friedmann-Robertson-Walker 
metric (1.3) 

ds'^ = dt^ - a^it) dx^ , so that dx = ■ (4.4) 

a[t) 

Using the definition of the redshift (1.9), we can turn the dependence on time 
into a dependence on the redshift: 



dz d /ao\ dQ , ttQ 
dt dt \ a / a 



^d = -^H. (4.5) 



The Hubble parameter is described by the Friedmann equation (1.14). In 
accordance with observations we assume a spatially fiat matter- and cosmo- 
logical constant-dominated universe with Q/^ + = 1 and k ~ Q/^/Qjn ~ 3. 
These parameters describe very well the cosmology of our universe up to 
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z = (9(1000), where the radiation content plays no dominant role. In this 
case, the Hubble parameter is given as 



Hiz) = Hoil + zf ^^]„(l + «:(l + z)-=^) (4.6) 
and the differential of the comoving distance becomes 

(l + ^)-3/2dz 

dx = , — . (4.7) 



aoHo^nm{l + n{l + z) 

Consequently, we obtain for the extragalactic contribution to the differ- 
ential neutrino flux per unit energy 

dJ f dN 

'-^e, / / =, (4.8) 



dE J U{Ey)./i 



3 



where y = 1 + z and z^ec is the time when the neutrinos decouple from the 
thermal plasma and start to propagate freely.^ The coefficient is given by 

VLdmQc 

A, 



1.1 X lO-^fcm^ssr)-! ( ''"\^ V f^^V' ■ 

1.3 X 1026 s/ VlSOGeV/ 



To obtain this result we have taken the gravitino density to equal the cold 
dark matter density, Vt^i2h'^ = QoAih"^ — 0.1, and the other constants as 
gc ^ 1.05 X 10"^ GeVcm-^ fi^ ~ 0.25 and h ^ 0.73. 

For a monochromatic spectrum, the redshift integral can be solved ana- 
lytically and results in [4] 



dE 1713/2 



1 + K 



2E ^ ^ 



"^3/2 



"V-)-%(,_ - I, (,10) 

V"^3/2/ V "^3/2 



where Q{x) denotes the Heaviside step function. This expression produces 
a characteristic 'triangular' contribution to the spectrum that terminates 
abruptly at the threshold. 



^In fact, the expression for the redshift integral is vahd only from the beginning of 
the matter dominated universe onwards. However, since the contributions from high z are 
negligible, we can actually perform the integration up to y = oo. 
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Halo Contribution In addition to tlie isotropic extragalactic signal there 
exists an anisotropic signal from gravitino decays in the Milky Way halo. 
Thanks to the small extension of the Milky Way compared to cosmological 
scales, we can work with distances in a fiat space-time. Since we are interested 
in the signal observed from Earth, it is convenient to use galactic coordinates, 
that have their origin in the Sun. In these coordinates spatial positions are 
specified by the distance s to the Sun, the galactic latitude b and the galactic 
longitude /. Analogous to the extragalactic case, the number of neutrinos 
from gravitino decay per unit energy and time in a volume element located 
at the coordinates / = (s, b, I) is given by 

^ = !^^§.^cosM.<iM/. (4.11) 

dEdt T3/2 dE ^ ' 

However, in this case the gravitino number density depends on the position 
in the galactic halo: 

(4.12) 

where we assume that the Milky Way dark matter halo consists only of 
gravitinos. The halo density distribution depends on the halo model but is 
spherically symmetric with respect to the galactic center in any case. A list 
of established halo models is given in Section 1.2. 

The distance r in the density profiles is given with respect to the galactic 
center, so we need to express it in terms of galactic coordinates. We find 

r(s, 6, /) = \J s"^ + Tq — 2 stq cos b cos / , (4-13) 

where ^ 8.5 kpc is the radius of the solar orbit around the galactic cen- 
ter [24]. Due to the far distance to the galactic center the observed signal is 
anisotropic. It is given by 

dJhaioib, I) dN^ 

^^ = A.„4^0^, (4.14) 
where the halo flux parameter is given by the line-of-sight integral 

^halo 



AhaioibJ) = / ghaio{r{s, b, l))ds . (4.15) 



In this expression Shaio is the extension of the dark matter halo in the cor- 
responding direction.^ Since the anisotropy is not very strong if the galactic 



^For the density profiles given in Section 1.2 the halo radius can safely taken to be 
infinite since the profiles are falling off at least proportional to for large radii. 
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center is excluded, we will use in the following analysis only an averaged 
signal. The averaged full-sky halo flux parameter is computed via^ 



1 



27r 7r/2 



Ahalo = ^ j I ^haioib, I) cosbdb. 





(4.16) 



and the averaged halo flux parameter excluding the galactic disk (galactic 



/ 2tt tt/2 



A 



excl.disk 



\0 7r/18 



is given 


by 










cos b db 


X 


27r y 


I 







• (4.17) 



It is also interesting to see how much the signal is increased in the direction 
of the galactic center. In order to perform an integration over a cone around 
the galactic center, it is convenient to change to new angular coordinates B 
and $ by the transformations 



cos b cos / — )• cos 9 , 

cos bdbdl — i- sin B dQ d^ . 



(4.18) 



The angle G describes the declination of the line of sight to the direction 
towards the galactic center and the angle $ gives the position in the rota- 
tionally symmetric direction around the connecting line between the Sun and 
the galactic center. The averaged halo flux parameter from a cone towards 
the galactic center with a half-cone opening angle of 5° can then be computed 
via 



A. 



GC 



57r/180 

j Ahaio{Q.^)^i^Qde 





X 



57r/180 

27r / sin 6 dQ 







(4.19) 

The numerical halo flux parameters obtained from different halo proflles 
are presented in Table 4.1. As expected from the different slopes in the inner 
parts of the Milky Way (see Figure 1.2), the results differ most for the cone 
directed towards the galactic center. By contrast, the results for the averaged 
halo flux parameters excluding the galactic disk agree within a few percent. 



■^Due to the symmetry of the halo profile we actually only integrate over one hemisphere 
in the galactic latitude. 
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Ahalo 




-^excl.dis 


k 


Agc 




NFW 


6.8 X 10- 


-8 


6.5 X 10 


-8 


4.2 X 10" 


-7 


NFWc 


7.1 X 10- 


-8 


6.4 X 10^ 


-8 


9.4 X 10" 


-7 


Moore 


6.9 X 10- 


-8 


6.4 X 10 


-8 


7.6 X 10- 


-7 


MoorCc 


7.2 X 10- 


-8 


6.4 X 10 


-8 


1.3 X 10- 


-6 


Kra 


5.8 X 10- 


-8 


5.6 X 10 


-8 


1.9 X 10- 


-7 


iso 


7.0 X 10- 


-8 


6.7 X 10- 


-8 


2.5 X 10" 


-7 



Table 4.1: Halo flux parameters for diflFerent halo models for the averaged 
full-sky signal, the signal excluding the galactic disk and the signal from the 
galactic center, respectively. The numerical values for the parameters are given 
in units of (cm^ s sr)-i(1.3 x 10^6 s/r3/2)(150 GeV/mg/s). 



Thus, in the following sections we will only present the results obtained 
using the Navarro-Frenk- White profile that is given by (cf. Section 1.2) 

QNFwir) = , (4.20) 

(r/rc) [1 + (r/rc)J 

with Tc = 20kpc and a normalization qq = 0.26GeVcm-^. The halo flux 
parameter for the averaged signal excluding the galactic disk is given, in this 
case, by 

.NFW 8/ 2 , 1 /1.3 X 1026s\ /150GeV\ 
^S... = 6.5xl0-«(cm2ssr)-M — . (4.21) 



T3/2 J V "^3/2 / 

This value will be used in most parts of the further analysis 



4.1.2 Neutrino Propagation and Neutrino Oscillations 

After being produced in gravitino decays, the neutrinos propagate over astro- 
physical distances until they can be observed on Earth. While propagating 
the neutrinos undergo flavor oscillations. Following the review in [53] we will 
thus shortly summarize the relevant consequences of neutrino oscillations. 

The neutrino gauge eigenstates Ue, and z/^ do not coincide with the mass 
eigenstates i/j, i = 1, 2, 3. Since the neutrino masses are negligible but dif- 
ferent from each other, there arise oscillations between the different neutrino 
flavors. The quantum mechanical treatment of neutrino oscillations leads to 
formulae for the probability that a neutrino with specific flavor is converted 
into a neutrino with another flavor. There are several general constraints for 
these probabilities: CPT-invariance implies 

Pii^i ^ u,) = P{vv ^ 1^1) . (4.22) 
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If we assume CP-invariance (and thus also T-invariance), we have 

The conversion probabihties for the different flavors are then given by [53] 
P(z/e ^ z/^J = S23 sin^ 26'i3S'23 + C23 sin^ 2612812 , 

P(z/e ^ Vr) = C23 Sin^ 2^13^23 + S^3 slu^ 2^12^12 , (4.24) 

P{v^ ^ Vr) = Ci3 sin^ 26'235'23 - S23C23 sin^ 26'i2S'i2 , 
while the survival probabilities are 

P(z/e Ue) = 1 - Sin^ 26'i3S'23 - 0^3 slu^ 26'i2S'i2 , 

P{u^ ^u^) = l-A 45^3(1 - 45^3)523 - 4 sin^ 2^12^12 , (4.25) 

P(z/^ U^) = 1 - 4Ci3C23(l - Ci3C23)S'23 " S23 siu^ 26'i25'i2 • 

In the above expressions we used the abbreviations 

Sij = sin 9ij , Cij = cos 9ij and Sij = sin I — . (4.26) 

\ AE J 

Using convenient units, we have for the last expression 

Am?.(eV2)L(km)\ ^ ^ 

(^-^^ E(GeV) )• ^'-''^ 

From the roots of this function we can infer the oscillation length for two- 
fiavor oscillations: 

AttE , E(GeV) 

= - 2-48 km ^ ; , J^^ . 4.28 

^ Ami Am2. (eV2) ^ ^ 

The experimental best-fit values for the neutrino mixing angles 612, O23 
and ^13, and the neutrino mass differences Am^^- = — mf are given in 
Table 4.2. From the value of ^23 we see that there is maximal mixing between 
muon and tau neutrinos. 

The cosmological dark matter density is spatially constant and the dark 
matter density in the galactic halo varies only on scales that are large com- 
pared to the oscillation lengths at energies in the GeV range. This leads to 
an averaging of the neutrino oscillations from the viewpoint of the observer. 
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sin^^ia sm2^23 sin^ ^13 Am22 (eV^) jAmfgl (eV^) 
0.304 0.50 0.01 7.65 x 10"^ 2.40 x 10"^ 

Table 4.2: Neutrino oscillation parameters. Figures taken from [54]. 

The huge path length combined with the finite energy resolution of detec- 
tors adds to this effect. Therefore, we are in the so-called long baseline limit, 
where 5*12 = 5*23 = 1/2. 

Neglecting the small mixing angle 6*13, the resulting oscillation probabili- 
ties are 

P{Ue ^ Ue) = 0.56 , 

P(z/e ^ U^) = P{Ue ^ = 0.22 , (4.29) 

Pij^t. ^ Jy^) = P{y^ ^ Vr) = P{Vr ^ T^r) = 0.39 . 

All other oscillation probabilities are determined by the relations (4.22) and 
(4.23). We note that even when the primary neutrino fiux is originally mainly 
composed of tau neutrinos, the flavor oscillations during the propagation 
produce comparable fluxes in all neutrino flavors. In particular, due to the 
maximal mixing between muon and tau flavor, the signals for the muon and 
tau neutrinos are identical after propagation. Thus, the initial assumption 
that the i?-parity breaking coupling is generated through a sneutrino VEV 
in the tau generation is not crucial. 

The fluxes for the different neutrino flavors and their extragalactic and 
halo contributions are shown in Figure 4.1. As mentioned before, we use the 
averaged halo flux excluding the galactic disk. Additionally, we employ an 
energy resolution of AE/E = 10%. Therefore, we convolute the neutrino 
flux with a Gaussian distribution 



^{E,E') = e (4.30) 

v27r a 



in the following way: 

dJi/ Gauss /" .T,/ 7-1 r-i/\ (^Jv 



dE 



I <l>{E,E')^dE'. (4.31) 



The standard deviation of the Gaussian distribution is chosen according to 
the energy resolution: a = AE/E. 

Note that even with this optimistic assumption for the energy resolution 
the lines from the decay into 7 z/,- and h Vj. become practically indistinguish- 
able from the continuum whereas the line from the decay into Z^v^ can be 
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1 10 100 

E (GeV) 

Figure 4.1: Total neutrino fluxes and the extragalactic and halo contributions 
for the different neutrino flavors after propagation using an energy resolution 
of 10 %. The gravitino mass and lifetime are chosen to be 7713/2 = 150 GeV and 
r3/2 = 1.3 X 10^6 s. 

resolved. Its position could allow a determination of the gravitino mass, even 
without determining the endpoint of the spectrum. 

In addition, we note that the extragalactic contribution to the neutrino 
fluxes is subdominant. Nevertheless, we will include this contribution in the 
numerical results. 

4.2 Neutrino Background 

The detection of a possible neutrino signal from gravitino decay is hindered 
by considerable neutrino backgrounds. Namely, in the energy range of 1- 
1000 GeV, there exist large background neutrino fluxes produced by interac- 
tions of cosmic rays with the Earth's atmosphere or with the solar corona, 
as well as neutrino fluxes from distant galactic sources. Let us briefly discuss 
these backgrounds separately. 

4.2.1 Atmospheric Neutrinos 

Collisions of energetic cosmic rays with the Earth's atmosphere initiate cas- 
cades that generate atmospheric neutrinos. When the primary cosmic rays 
hit the nuclei of air in the upper part of the atmosphere, mostly pious and 
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some kaons are produced. The charged pions decay dominantly into muons 
and muon neutrinos via 

tt"^ li^i'^ and 7r~ IJ'~i^n , 

while the decay into electrons and electron neutrinos is suppressed by ml/rnj^. 
Before decaying, the muons produced in the pion decays travel an average 
distance 

4-cr,7.-lkm^^, 

where r^u is the muon lifetime and 7^ = E'^/m^i is the relativistic dilatation 
factor. If the muons decay via 

before they hit the Earth's surface, one obtains a flux ratio of 2:1 for muon 
and electron neutrinos and antineutrinos. However, for energies above a few 
GeV the muons hit the surface before decaying. Therefore, at higher energies 
the ratio {p^ + v^jiVf, + is larger than 2. 

In order to determine the fluxes of atmospheric neutrinos one has to 
consider several effects. First of all, the spectrum of cosmic rays and their 
chemical composition has to be known. For energies above a few GcV the 
cosmic ray flux has an energy dependence d(j)/dE oc E~^''^ and most of its 
constituent particles are protons. At low energies (< 10 GeV) the cosmic 
ray flux is modulated by the solar wind and is therefore dependent on the 
eleven-year-cycle of solar activity. 

Another important effect on lower energy cosmic rays comes from the 
geomagnetic field. Since the cosmic ray particles travel along the field lines 
at the poles and perpendicular to the field lines at the equator, this effect 
is dependent on the position on Earth. So there exists a position-dependent 
geomagnetic cutoff in the magnetic rigidity TZ — p/{Ze), where p is the 
momentum and Ze the charge of the cosmic rays. Particles with a magnetic 
rigidity below the cutoff do not reach the atmosphere. 

These effects and the details of the interactions and decays in the atmo- 
sphere can be incorporated into a Monte Carlo procedure, making it possible 
to calculate the atmospheric neutrino fluxes with high accuracy. The largest 
uncertainties stem from the initial cosmic ray spectrum and the hadronic 
interactions, leading to a ~ 20 % uncertainty in the absolute neutrino fluxes. 
However, the ratios of the fluxes of the different neutrino flavors are accurate 
to better than 5 %. 

At energies above 1 GeV the fluxes of the atmospheric neutrinos obey an 
approximate power law: 
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for electron and muon neutrinos, respectively. At lower energies the slope is 
shallower due to geomagnetic effects and solar modulation. 

In this work we use the fluxes of atmospheric electron and muon neutri- 
nos computed by Battistoni et al. with the Monte Carlo package FLUKA 
(Fluktuierende Kaskade) assuming massless neutrinos [55]. For the energy 
range 0.1-112 GeV we use the fluxes for the Kamioka site including the ge- 
omagnetic cutoff [56]. To be conservative, we use the larger fluxes from the 
minimum solar activity sample. For the energy range 25-10000 GeV only the 
muon neutrino flux is available [57]. 

The main tau neutrino background stems from atmospheric muon neu- 
trinos that oscillate into tau neutrinos on the way to the detector. From 
equation (4.28) we see that the oscillation length on Earth is too short for 
a significant conversion of electron neutrinos. Therefore, we apply two-flavor 
neutrino oscillations between muon and tau neutrinos to the initial fluxes. 
Since the mixing between muon and tau neutrinos is maximal (sin 26*23 — 1), 
the conversion probability is given by 



In this expression E is the neutrino energy and L their propagation length 
after the production in the atmosphere. The latter depends on the zenith 
angle and is given by 



with ~ 6.4 X 10^ m being the mean Earth radius and h 15 km being 
the mean altitude at which the muon neutrinos are produced [58]. Moreover, 
the relevant neutrino parameters for atmospheric oscillations are given in 
Table 4.2. 

4.2.2 Further Neutrino Background Sources 

In addition to the flux of tau neutrinos originating from the conversion of 
muon neutrinos, there exists an intrinsic contribution from the decay of 
charmed particles that are produced in the atmosphere in cosmic ray in- 
teractions. This flux is suppressed by a factor of about 10^ compared to the 
flux of electron and muon neutrinos from pion and kaon decay. 

This intrinsic contribution has been computed by Pasquali and Reno [59] 
in a PYTHIA based calculation and can be parameterized as 




(4.33) 



L = ^/{R^coseY + 2R^h + h^ - cos 6 , 



(4.34) 



'dE 



(4.35) 
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Flavor 


No 


7 


A 


^0 (GeV) 


7' 




Ue + I'e 

+ 


7.4 X 10-6 
1.3 X 10"^ 


2.03 
1.98 


8.5 X 10-*^ 
8.5 X 10-*^ 


1.2 X 10^ 
3.0 X lO'' 


2.33 
2.38 


5.0 X 10"^ 

5.1 X 10-3 



Table 4.3: Values for the parameterization of the corona electron and muon 
neutrino flux. 



where x = log^Q E (GeV), A = 6.69, B = 1.05, C = 0.150 and D = -0.00820. 
The next-to-leading-order perturbative QCD calculation also presented in the 
paper gives lower fluxes at energies below several TeV and is therefore less 
conservative. 

Analogous to the production of neutrinos in the Earth's atmosphere, neu- 
trinos are produced in the solar corona by cosmic ray collisions. This neutrino 
flux has been studied by Ingelman and Thunman in [60], who found that the 
flux of electron and muon neutrinos integrated over the solar disk can be 
described by the following parameterization: 

GeV- cm- s- „,'::1!S-.' _ „ (4-36) 



dE 1 KiEjGeY))- 

1+AE (GeV) ' 



E>E, 







It is valid for 10^ GeV < E < 10^ GeV and the numerical values of the 
parameters are given in Table 4.3 for both neutrino flavors. 

The electron and muon neutrinos and antineutrinos produced in the solar 
corona oscillate during their propagation to Earth. Considering the large 
distance traveled (compared to the oscillation lengths), the conversion and 
survival probabilities can be averaged, and the fluxes on Earth in the different 
flavors can be obtained from equation (4.29). 

Last, the fluxes of tau neutrinos that originate from galactic sources are 
discussed by Athar, Lee and Lin in [61]. For the tau neutrino flux from the 
galactic plane in the presence of neutrino oscillations they flnd the parame- 
terization 

% = 9 X 10-6 GeV-i cm-2 s'' si'' {E (GeV))"''^' , (4.37) 
dE 

that is valid in the energy range 1 GeV < E < 10^ GeV. 



4.3 Reduction of the Neutrino Background 

The averaged full-sky signal for the neutrinos from gravitino decay is shown 
in Figure 4.2 together with the results for the atmospheric background from 
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Figure 4.2: Averaged fuh-sky neutrino fluxes expected at the Super-K site. 
The gravitino parameters and the energy resolution are chosen as in Figure 4.1. 



FLUKA. The signal lies several orders of magnitude below the expected atmo- 
spheric background for all flavors. Therefore, we find that the interpretation 
of the EGRET and HEAT anomalies in terms of gravitino decay is compati- 
ble with neutrino flux measurements, as it does not lead to an overproduction 
of neutrinos. 

Going beyond this consistency check, we will examine in the following 
section the possibility of detecting this extraordinary flux contribution in 
neutrino experiments. Due to the low signal-to-background ratio, the signal 
cannot be detected directly. It will therefore be necessary to find strategies for 
effectively reducing the background in order to have any chance of detecting 
the signal. As is apparent from Figure 4.2, the tau neutrino channel appears to 
be the most promising of the three flavors, since it has the lowest background. 

In general, the neutrino spectrum from gravitino decay has some very 
characteristic features that could allow to distinguish it from the featureless 
backgrounds, but the question is whether neutrino detectors are able to reach 
sufficient sensitivity to resolve these features. 



4.3.1 Electron and Muon Neutrinos 

For the electron and muon neutrinos, that are more easily detected in neu- 
trino observatories than tau neutrinos, the signal-to-background ratio is very 
small (~ 10~^-10~^) for gravitino lifetimes that are not already excluded by 
bounds from the gamma-ray [5, 24] or antimatter signals [6, 7]. As can be 
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seen from Figure 4.2, even the peak of the spectrum is three orders of magni- 
tude below the background. This would make distinguishing an extraordinary 
signal from the background appear extremely difficult. 

Unfortunately, we could not find a suitable strategy to sufficiently reduce 
this background, for example by exploiting directionality. In general, the at- 
mospheric neutrino baseline (L < 2 R^) is too short for a complete conversion 
of muon neutrinos into another flavor at energies on the order of 50-100 GeV 
(cf. equation (4.28)). Therefore, a detection of the signal without having 
prior knowledge of the position of the peak in the neutrino spectrum from 
gravitino decay seems hopeless. If information on the spectrum is available, 
for instance from the detection of a monochromatic gamma-ray line in the 
Fermi Gamma-ray Space Telescope (FGST, formerly known as the Gamma- 
ray Large Area Space Telescope (GLAST) [62]), one could envisage strategies 
to disentangle the neutrino signal from the background. Anyway, that would 
probably require a much better knowledge of the atmospheric neutrino flux 
at the relevant energies and a better energy resolution of neutrino detectors 
than presently available. 



4.3.2 Tau Neutrinos 

For the tau neutrinos the signal-to-background ratio is more promising since 
it lies above ~ 10^^ at the peak energy (cf. Figure 4.2). Moreover, most of 
the tau neutrino background from atmospheric oscillations can be effectively 
reduced by exploiting directionality. Since we are interested in energies of 
(9(10-100) GeV, we are in the region where the oscillation length for two- 
flavor oscillations between muon and tau neutrinos is larger than the Earth 
diameter. Thus, the conversion probability for down-going muon neutrinos 
(i.e. neutrinos from above the horizon) is approximately given by 

(..^^^Mgp)^. ,,33, 

We have seen earlier that the atmospheric muon neutrino flux has an energy 
dependence oc in the range above 1 GeV. Therefore, we expect for the 
down-going tau neutrinos a flux with an energy dependence 

oc E-^. (4.39) 



dE 



Due to the steeper slope, the previously subdominant tau neutrino back- 
ground sources become important at higher energies. Anyway, these back- 
grounds can also be effectively reduced by exploiting directionality. In fact. 
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Figure 4.3: Down-going tau neutrino fluxes expected at the Super-K site 
with (left) and without (right) the contribution from the galactic disk and the 
solar corona. The gravitino parameters and the energy resolution are chosen 
as in Figure 4.1. 

since the solar corona neutrinos and the galactic neutrinos come from the 
direction of the Sun or from the galactic plane, these directions could sim- 
ply be excluded from the search to reduce the background. This is why we 
excluded the galactic disk in the calculation of the halo contribution to the 
neutrino flux. 

Thus — regarding only down-going tau neutrinos — the background can be 
reduced by several orders of magnitude. In Figure 4.3 we show the fluxes for 
down-going tau neutrinos at the Super-Kamiokande site. We see that in this 
case the signal can overcome the simulated background from FLUKA, even 
without cutting away the Sun or the galactic plane. 



4.4 Detection Prospects 

In this section we discuss the prospects for the detection of the extraordinary 
signal in the tau neutrino flavor in present and future neutrino experiments. 
For that purpose, we first briefly describe the detection technique of water 
Cerenkov detectors following the reviews in [53, 63]. 

4.4.1 Water Cerenkov Detectors 

Water Cerenkov detectors like Super-Kamiokande or IceCube consist of pure 
water or ice and a large number of photomultiplier tubes (PMTs) that detect 



63 



Chapter 4. Neutrinos as a Probe of Gravitino Dark Matter 



0.6 




1 10 100 

E (GeV) 

Figure 4.4: Total cross sections of CC interactions for Vr (solid) and z^r 
(dashed) calculated in [64] using the program NEUT. 

in real time the Cerenkov light emitted by relativistic charged particles. 

Charged leptons can be generated by neutrinos, that cross the detector, 
via charged current (CC) scatterings off nuclei in the water {viN ^ I N') and 
in the case of electron neutrinos also by elastic scatterings off electrons. Fol- 
lowing these interactions, Cerenkov light is emitted if the charged particle's 
velocity is larger than the speed of light in water. 

This technique works well for electron and muon neutrinos, which gen- 
erate an electron and a muon in the CC interaction, respectively. Muons 
produce a sharp Cerenkov ring, while electrons produce a diffuse ring. Thus, 
the neutrino flavor can be identified using the Cerenkov light signature. On 
the other hand, it is not possible to distinguish particles from antiparticles. 

Tau leptons, generated from tau neutrinos, decay very fast via hadronic 
channels and therefore do not produce a clear Cerenkov light signature. The 
cross section for CC interactions of tau neutrinos and antineutrinos with 
nuclei has been calculated with the neutrino interaction simulation program 
NEUT in [64]. The result is presented in Figure 4.4. Clearly visible is the 
threshold for r lepton production in the reaction t/^ n — )■ r p at E^^ > rur + 
m2/2m„ ~ 3.5 GeV. 

By the measurement of the Cerenkov light with PMTs, it is possible to 
reconstruct the energy and the direction of the charged lepton. However, 
this is not sufficient to reconstruct the neutrino energy without knowing the 
direction of the incoming neutrino, which differs from the lepton direction 
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by A6' ~ 30°^/l/E,y (GeV) [65]. Therefore, the neutrino events are usually 
classified according to the signal topology. 

Electron or muon events for which the particle track starts and ends inside 
the detector are called fully contained events. In this case, the lepton energy 
can be reconstructed. This event type is further divided into sub-GeV and 
multi-GeV events, depending on whether the lepton energy is below or above 
~ 1.4 GeV. On average, sub-GeV events originate from neutrinos of several 
hundred MeV, and multi-GeV events arise from neutrinos of few GeV energy. 

If the lepton is produced inside the detector but escapes from it, the 
event is called partially contained. In this case, only the visible energy of 
the lepton can be reconstructed. Since these events originate from neutrinos 
with energies comparable to those of fully contained events, they are usually 
grouped together. 

When a muon is produced in the rock below the detector and its track 
ends inside the detector, one speaks of up-going stopping muon events. The 
typical energy of parent neutrinos is ~ 10 GeV. However, this event topology 
does not occur for electrons, as they produce electromagnetic showers in the 
surrounding material before reaching the detector. We note that down-going 
muons are usually not taken into account, since they cannot be distinguished 
from the background of cosmic ray muons. 

If a muon track crosses the whole detector from below, the event is called 
up-going through-going muon. These events originate from neutrinos with 
energies of 10-1000 GeV. 

The contribution of neutrinos to the different event topologies described 
above can be simulated and is presented in Figure 4.5 as a function of the 
neutrino energy. From this discussion we see that it is only possible to re- 
construct the neutrino energy on a statistical basis using the event topology 
and the constraint that the neutrino energy must be larger than the visible 
energy of the charged lepton. 

Given this fact, the employed energy resolution of 10 % for the neutrino 
fiux in all figures is very optimistic and way beyond feasibility in present 
neutrino detectors. 

4.4.2 Observability in Super-Kamiokande 

Super-Kamiokande is a cylindrical 50 kton water Cerenkov detector located 
in the Kamioka mine in central Japan. The fiducial mass of the detector is 
22.5 kton (i.e. the effective detector mass that can be used for analyses). 

The Super-Kamiokande collaboration has developed a statistical method 
to discriminate tau neutrinos from the background of other flavors [66]. Using 
two different strategies, namely a likelihood analysis and a neural network. 
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Figure 4.5: Contributions of neutrinos to distinct event topologies as a func- 
tion of the neutrino energy. Figure taken from [63]. 

they find an efficiency of 43.1 % and 39.0 %, respectively, to identify tau neu- 
trinos correctly. However, they still misidentify 3.8 % and 3.4 %, respectively, 
of the electron and muon background neutrinos as tau neutrinos. 

Due to the large number of atmospheric electron and muon neutrino 
events, the sample of tau neutrinos is thus dominated by misidentified neu- 
trinos. The actual tau neutrino events can therefore only be extracted on a 
statistical basis, using Monte Carlo methods. In the end, the data is found to 
be consistent with the expected atmospheric tau neutrino flux in the presence 
of neutrino oscillations: The full-sky atmospheric tau neutrino signal results 
in fact in 78 events from CC interactions in the fiducial volume during the 
Super-K I period and 43 events during the Super-K II period [64, 66]. 

However, the above mentioned analysis does not exploit the information 
about the spectral shape of the signal, apart from setting a threshold for r 
lepton production. So this kind of data analysis could certainly be improved 
in order to search for a signal with a peak above the continuum, as in our 
case. 

Apart from the experimental difficulties, it is worthwhile to examine the 
theoretically expected signal in the tau channel. Figure 4.6 shows the ex- 
pected number of tau neutrino and antineutrino events per century of obser- 
vation at Super-Kamiokande within a zenith angle integrated from cos 6' to 1. 
If only down-going neutrinos are selected, the signal from the gravitino decay 
lies above the atmospheric background, especially for higher energies. How- 
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Figure 4.6: Number of events per century of observation at Super-K due 
to CC interactions of Vr and Vr from the atmosphere (dot-dashed) and the 
gravitino decay (sohd), integrating the signal from the zenith direction to the 
angle 9. The gravitino mass and lifetime are chosen as in Figure 4.1. 

ever, the fluxes are extremely low and result only in a few events per century, 
making it practically impossible to discriminate them from the other flavors 
using statistical methods. 

One detector specifically optimized for measuring tau neutrinos above 
17 GeV event by event is OPERA, which is already active in Gran Sasso and 
will measure tau neutrino appearance in a muon neutrino beam produced 
at CERN. Unfortunately, the detector's effective mass is more than a fac- 
tor 10 smaller than that of Super-K and thus — even neglecting the issue of 
directionality — would be able to observe only one event in more than 1000 
years. 

We therefore conclude that present detectors are not able to detect the 
signal, either because they do not have sufficient efficiency for identifying tau 
neutrinos, or because they are too small for the low intensity of our signal, 
or both. 
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4.4.3 Future Detectors 

We will discuss briefly in this section the prospects of detection for the future. 
Hyper-Kamiokande 

The prospects for Hyper-Kamiokande can be easily obtained by considering 
that its mass is planned to be a factor of 10 (for the O.SMton project) to 
20 (for a 1 Mton case) larger than Super-K. Assuming that the rest of the 
detector characteristics are unchanged, we expect it to have at most 0(20-40) 
events from our signal per century from the upper hemisphere.^ 

This number of events might be still too small to allow for statistical 
analysis. On the other hand, we expect most of the events to appear within 
the peak region or near the threshold. Therefore, using an appropriate energy 
binning — especially if the signal has been detected already in the gamma-ray 
channel — could even allow to collect a significant number of events above the 
background in a specific energy bin on shorter timescales. 

Nevertheless, it is clear that a sufficiently good energy resolution is a key 
requirement for extracting the events from the line, and it remains uncertain 
how and if the statistical tau flavor discrimination analysis can be applied to 
a sample of such few events. 

IceCube and km^ Detectors 

Detectors of km^ dimensions have in principle sufficient size to collect enough 
events to detect the signal from gravitino decay within a reasonable time 
span. Even considering that IceCube is actually looking downwards and not 
at the upper hemisphere, from the horizontal direction and the CC cross sec- 
tion for tau neutrinos we estimate (9(100) events per year for the completed 
experiment. Of course, the effective area depends on the neutrino energy: 
Taking the effective area, given in [67] for up-going events, for our signal from 
the opposite direction and assuming most of the signal is above 100 GeV, we 
have instead C(10) events per year. 

In general, it would be desirable to decrease the energy threshold of 
IceCube below 100 GeV in order to cover the energy range favored by the 
EGRET and HEAT anomalies. The combination of IceCube with the Antarc- 
tic Muon And Neutrino Detector Array (AMANDA) already allows to lower 
the threshold to 30 GeV. Additionally, plans are being considered for adding 
another, denser subdetector at a deeper location to improve the sensitivity 

^This number could be increased if the fiducial volume of Hyper-Kamiokande is larger 
than ~ 1/2 of the total volume, as it is in Super-Kamiokande. 
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to dark matter annihilations [67]. Such a configuration could probably also 
be useful for investigating the present scenario and, more generally, other 
decaying dark matter candidates. 

As we have stated before, in the case of water Cerenkov detectors, the dis- 
crimination of tau neutrinos from other neutrino fiavors is generally difficult, 
and for IceCube strategies for tau fiavor identification have been proposed 
only for neutrinos well above TeV energies [68]. It could therefore be more 
favorable to improve the energy resolution and exploit the muon neutrino 
final state instead. 



4.5 Constraints on the Gravitino Parameters 

In this final section we want to present constraints on the gravitino lifetime 
as a function of its mass from the non-observation of the neutrino signal from 
gravitino dark matter decay. To achieve this, we compare the neutrino fiux 
from gravitino decay to the atmospheric background fiux. 

In particular, we only take the halo contribution of the photon and 
lines discussed in Section 3.2 to approximate the signal. 



excl.disk -"-^^v / "ry ^ 2 J 



dJ: 



(4.40) 



dE — ^ " {E^-ElY + ElVl 

and consider either the signal in all neutrino fiavors or only the tau fia- 
vor. According to the unitarity of the conversion probabilities for neutrino 
oscillations, the signal in all neutrino fiavors is given by the expressions 
above, whereas for the tau signal we have to apply the survival probabil- 
ity P{i'r I^t)- 

To allow for the poor energy resolution of neutrino detectors, we adopt a 
nominal energy resolution of a = 0.3 in log^g (^eV) around the peak posi- 
tion [65]. Thus, we convolute the neutrino fiuxes with a Gaussian distribution 
in the logarithm of the neutrino energy: 



2 



^{E,E') = ^^---e ^ ^ (4.41) 

y l-KG In 10 -c/ 



where the altered coefficients take care of the correct normalization and the 
convolution is computed according to 

oo 

dJy,Qa,uss f 1/7-1 7-i/\ dJi, 



dE 



I ^E,E')^dE'. (4.42) 
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Figure 4.7: Region of the gravitino lifetime where the hne signals from the 
two-body decays into 7 v-r- and Z^Ut overcome the atmospheric background 
either for all neutrino flavors or for down-going tau neutrinos only. We see 
that considering the tau flavor and exploiting directionality can improve the 
constraints by three orders of magnitude for large gravitino masses. We use 
here the branching ratios shown in Figure 3.1 for the two lines. Note that in the 
region above the threshold the continuum spectrum from its fragmentation 
may be used to close the gap, but we do not consider this possibility here. 



The actual lower bound on the gravitino lifetime is determined by the 
criterion that the flux at the peak position equals the simulated background 
either for all flavors or for down-going tau neutrinos only. The result is pre- 
sented in Figure 4.7. For the flux in all neutrino flavors, it is similar to that 
presented in [65], except that we are weighting the channels with the gravi- 
tino branching ratios and that we have only one neutrino produced in the 
lines instead of two. 

Note that since the signal from gravitino decay is proportional to I/T3/2, 
requiring the peak flux to be larger than the background by a particular 
factor only rescales the exclusion curves by the inverse of this factor. 

We clearly see that again the tau neutrino channel in the down-going 
direction allows to constrain the gravitino lifetime a few orders of magnitude 
better than the muon or electron neutrino channels. On the other hand, 
similar plots for the gamma-ray channel are even more sensitive and give 
bounds on the order of 10^'' s for gravitino masses below the and Z*^ 
thresholds [24]. 

For masses above (9(100) GeV, the tau neutrino channel even starts to 
compete in sensitivity with the photon channel. Admittedly, we must point 
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out that the presented bounds are only hypothetical since we neglect the 
difficulties connected with the measurement of such a low flux in neutrino 
detectors and the identification of the neutrino fiavor. 
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The scenario of gravitino dark matter with shghtly broken i?-parity is the- 
oretically well motivated as it leads to a consistent thermal history of the 
universe including both baryogenesis via thermal leptogenesis and standard 
Big Bang nucleosynthesis. 

Although the above scenario was originally devised to reconcile the ther- 
mal history of the universe with supersymmetric dark matter, it has been 
pointed out in the recent literature that the decay of gravitino dark matter 
with a lifetime of ~ 10^^ s and a mass of ~ 150 GeV into massive gauge 
bosons may account for the anomalies observed in the diffuse extragalactic 
gamma-ray spectrum as measured by EGRET, as well as in the positron 
fraction as measured by HEAT.^ 

In this thesis we examined the neutrino signal from the decay of unstable 
gravitino dark matter in a scenario with bilinear i?-parity breaking. Moti- 
vated by the possible explanation of anomalies in other cosmic ray channels, 
we recomputed the gravitino decay widths, taking into account also the Feyn- 
man diagrams with the non-abelian 4-vertex for the and Z° channels 
and the decay into the lightest Higgs boson. Using the resulting branching 
ratios for the gravitino decay channels, we determined the neutrino signal on 
Earth originating from gravitino dark matter decay. 

As a consistency check of the scenario, we employed the same gravitino 
parameters as those that were used to account for the anomalies in other 
cosmic ray channels and found that the neutrino signal is compatible with 
results from present neutrino experiments. In addition, we have examined 
the detectability of this extraordinary component of the neutrino flux to find 
an independent way to test this scenario. 

^The existence of an excess in the positron fraction seems to be supported by prelim- 
inary results of the PAMELA (a Payload for Antimatter Matter Exploration and Light- 
nuclei Astrophysics) satellite mission [69]. 
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While the signal in the neutrino spectrum with two or more distinct peaks, 
resulting from two-body gravitino decays into gauge/Higgs boson and neu- 
trino, is very characteristic, it will be challenging to detect these features in 
neutrino experiments. On one side, present neutrino detectors do not achieve 
a sufficient energy resolution to resolve peaks in the spectrum, and on the 
other side, the event rate is expected to be so small that the background of 
atmospheric neutrinos overwhelms the signal in all flavors. 

The most promising signal-to-background ratio is found for tau neutrinos, 
especially when considering only the flux from above the horizon, since there 
the atmospheric tau neutrino flux is substantially reduced. However, tau neu- 
trinos are difficult to identify in water Cerenkov detectors and probably only 
an event-by-event identification could allow detecting the signal with such 
extremely limited statistics. Therefore, we conclude that it is currently not 
possible to detect this contribution due to present technological limitations. 

The ideal detector for testing the present scenario would be one of mega- 
ton mass with the ability to identify and measure tau neutrinos event by 
event. Should such a detector ever become available, it could be worthwhile 
to look for this extraordinary component of the neutrino flux by employing 
strategies for background reduction such as the ones discussed here, espe- 
cially if the anomalous signatures in the positron fraction and the diffuse 
extragalactic gamma-ray spectrum are conflrmed by PAMELA and FGST, 
respectively. 

In principle, the non-observation of the neutrino signal can be used to infer 
bounds on the gravitino signal, but we found that these bounds presumably 
cannot compete with the bounds from other decay channels, like the decay 
into photons, positrons or antiprotons. 

On the other hand, the detection of a signal in neutrinos, that is com- 
patible with signals in the other indirect detection channels, would in fact 
bring signiflcant support to the scenario of decaying dark matter, possibly 
consisting of gravitinos that are unstable due to bilinear i?-parity breaking. 
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Appendix A 

Units and Physical Constants 



Here we will shortly summarize the values of all the physical and astrophysical 
constants that appear in the thesis. At times, we work in units where the 
reduced Planck constant, the speed of light and the Boltzmann constant 
obey h = c = k = l.ln this case, we have the conversion factors 



1 eV = 1.160 4505(20) x 10^ K , 
leV = 1.602176487(40) x 10~^^ J, 
1 GeV = 1.782 661 758(44) x 10"^^ kg , 
IGeV"^ = 1.973 269 631(49) x lO^^^m, 
1 GeV"^ = 6.582 118 99(16) x 10"^^ s . 



All figures are taken from the The Review of Particle Physics [17], if not 
marked otherwise. Numbers in parentheses represent the one standard devi- 
ation uncertainty in the last digits. 



Quantity 



Symbol 



Value 



muon mean life 
tau mass 
tau mean life 
boson mass 

boson mass 

boson decay width 
neutron mass 
neutron mean life 
proton mass 
gravitational constant 



muon mass 



electron mass 



777.7- 



m, 



510.998 910(13) keV 
105.658 3668(38) MeV 
2.197019(21) X 10-S 
1.776 84(17) GeV 




290.6(1.0) X 10-15 s 
80.398(25) GeV 



91.1876(21) GeV 
2.4952(23) GeV 



n 



939.565 360(81) MeV 
885.7(8) s 




938.272 013(23) MeV 
6.708 81(67) X lO'^^ GeV'^ 
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Quantity 



Symbol Value 



parsec 


pc 


3.085 6776 x 10^^ m 




Solar distance from GC 


^0 


8.0(5) kpc 




Earth mean equatorial radius 




o.oloiol X iU m 




local halo density 




0.1-0.7 GeVcm-3 




present-day CMB temperature 


Tn 
u 


2.725(1) K 




normalized Hubble constant 


h 


0.73(3) 




p^itip^^l dprmitv 


n 
f c 


1 053 68fll) X 10^5 h'^ 


CtpV cyt\ 


matter density 


Til 


0.128(8) 




baryon density 




0.0223(7) 




dark matter density 




0.105(8) 




dark energy density 




0.73(3) 




total energy density 


^tot 


1.011(12) 




baryon-to-photon ratio 


V 


6.12(19) X 10-1° 




solar neutrino mixing angle 


sin^ 6*12 


0.304(+22) 


[54] 


atmospheric mixing angle 


sin^ 6*23 




[54] 


third neutrino mixing angle 


sin^ 6*13 


0.010(1J?) 


[54] 


solar neutrino mass difference 


Am?2 


7.65(1^;^) X 10-5 


[54] 


atmospheric mass difference 




2.40(+}?) X 10-3 


[54] 
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Appendix B 

Notation, Conventions and 
Formulae 



9i^>y = 9 



(B.l) 



Four- Vectors and Tensors Lorentz indices are depicted by small Greek 
letters, e.g. // = 0, 1, 2, 3 . The metric of Minkowski space is chosen to be 

/+1 \ 
0-100 

0-10 

\ -1/ 

Lorentz indices of four- vectors and higher rank tensors are raised and lowered 
using the space-time metric: 

a^^g^^^a" and b'' ^ g''" K . (B.2) 

We use in all cases a sum convention for Lorentz indices: 



//=0 



(B.3) 



Antisymmetric Tensor The totally antisymmetric tensor £^^'^P'^ is chosen 
such that 

£°i23 = +1 _ (B.4) 
This implies for the totally antisymmetric tensor with lowered indices, 



^ixvpa — gfxa g^p gpf gaS ^ 



(B.5) 



that £oi23 = -1 ■ 

The totally antisymmetric tensor in three dimensions Eijk is chosen such that 



£l23 = +1 



(B.6) 
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Pauli Matrices The Pauli sigma matrices read 
and satisfy the identity 

cr-cTj = 5ij + iEijkCTk ■ (B.8) 

Gamma Matrices The Dirac gamma matrices form a ChfFord algebra 

{7^7^ = 2^?'^^ (B.9) 

In most cases it is convenient to suppress the spinor indices of the gamma 
matrices, so we do not write them here. Using the above algebra and the 
metric (B.l) it can be easily seen that 

(7°)' = 1 and (7^)' = -l. (B.IO) 

Hermitian conjugation of the gamma matrices yields 

(7'^)^ = 7V7°- (B.ll) 

The Lorentz indices of gamma matrices are raised and lowered by the metric, 
just as in the case of four-vectors. That implies 

70 = 7° and n = -i\ (B.12) 

Thus, we see that raising and lowering of Lorentz indices is equivalent to 
hermitian conjugation. In addition to 7°, 7^, 7^, 7^ , a fifth gamma matrix 
can be defined 

7' = ilWlW = W7^V7^ (B.13) 

Using the algebra (B.9) it can be easily shown that 7^ has the following 
properties: 

{7', 7^}=0, 

(7')' = 1, (B.14) 
(7^)^ = 7^ 

The matrix 7^ is used to define the chirality projection operators 

Pl = ^(1-7') and ^ 1 (l + 7^) . (B.15) 
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These operators possess the projector properties 

P| = Pl, PI^Pr and P^Pr^PrPl^O. (B.16) 

A very common object in theories of particles with half-integral spin is 
the contraction of gamma matrices with the particle momentum. Therefore, 
we use the Feynman slash notation 

^ = 7%. (B.17) 

In order to obtain Lorentz scalars, in the formulation of theories of particles 
with half-integral spin there occur products of adjoint spinors and spinors. 
The adjoint spinor and adjoint vector-spinor are defined as 

V; = VV and i^^ = ^ll', (B.18) 

respectively. 

Majorana Spinors Majorana particles are invariant under charge conju- 
gation. Thus, the Majorana condition reads 

X^X'' = Cx^- (B.19) 
Here we use the charge-conjugation matrix C that has the properties 

= c-\ 

C^ = -C, (B.20) 

with 

r, -! 1 forP =1,7^,7' 

'^'^ \ -1 forP =Y:(t'"' ' ^ ^ ^ 

In the above expression we use the definition a^"^ = j"^] , /i < u . 

Dirac Spinors The free Dirac equation 

{i$ - m) ijjix) = (B.22) 
has four linearly independent plane wave solutions: 

iP{x) = u\p) e-^^"^ and ^(x) = v\p) e'^'^ s = ±^ , (B.23) 
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where the spinors u and v have to obey the following constraints: 

{p-m)u''{p) = and + m) v'ip) = . (B.24) 
The normalization of the Dirac spinors is given by [70] 

u' (p) u'' {p) = 2m 5''' and v' (p) v'' (p) = -2 m 5''' , (B.25) 
and the polarization sums for spin- 1/2 fermions are given by [70] 

^M^(p) u\p) =f + m, 

s 

Y,v"{p) v'{p) 



(B.26) 



— m . 



where the sum is performed over the fermion polarizations s = ±| . 

Polarization Vectors Free massless spin-1 particles in the Lorentz gauge 
{d^Afj^ = 0) and free massive spin-1 particles obey the Proca equation [27] 

[d^d^ + M^)Af, = 0. (B.27) 

This equation has the solutions 

A^ = e^{p)e-'P-^, (B.28) 

where is a polarization vector. The condition d^A^^ = demands that 

P^e^ = , (B.29) 

reducing the number of independent polarization vectors to three. 

For massless spin-1 particles we have the freedom to make the additional 
gauge transformation 

A^-^ A'^ = A^ + df'A with d^df,A = 0. (B.30) 

The condition for A is required by the Lorentz gauge condition. This addi- 
tional gauge freedom reduces the number of independent polarization vectors 
for massless spin-1 particles to two. 

The polarization vectors are normalized according to [40] 

ej*e^> = -5^^', (B.31) 
while the polarization sums are given by [70] 
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where the sum is performed over the two polarization states of the massless 
spin-1 particle A = ±1, and by 

Y.4*(P)'tip) = -(9..-^), (B.33) 

where the sum is performed over the three polarization states of the massive 
spin-1 particle A = ±1, 0. 

Identities for Gamma Matrices In the calculation of squared ampli- 
tudes we encounter traces of gamma matrices. For these the following iden- 
tities can be derived from the algebra (B.9): 



Tr(l) 


= 4, 


Tr (odd number of 7s) 


= 0, 


Tr (7^^) 


= 4^?^^ 


Tr (7^7''7''7") 


= 4 {g^^r - 


Tr (7^) 


= 0, 


Tr (7^7''7^) 


= 0, 


Tr (7''7^7''7"7^) 





g^'g"" +g^"g''n, (B.34) 



We also encounter contractions of gamma matrices. Using the algebra (B.9), 
we obtain the following identities: 

T^'Tm = 4 , 

7^''7. = -2 7^ . . 

7^V7m = 4^7^'', ^ ■ ' 

Any complex 4x4 matrix M can be expanded in terms of a basis of 
gamma matrices, 

^ = {l,7^c^'^^7^^7'}, (B.36) 
according to the relation [71] 

M = ^ ^ Tr [Mr) Ti , (B.37) 



where Fj is the dual basis. The dual basis can be determined using the or- 
thogonality relation 

Tr (FiF^) = 45/ (B.38) 
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and corresponds to the hermitian conjugate of the basis. Thus we have 

r, = {1, 7^, a^,, 7^7^, 7^} . (B.39) 

We want to apply this expansion identity to a combination of gamma 
matrices that appears in the couphng of the gravitino to the gauge super- 
multiplet: 

r [7^ Y] = I {Tr (7^^ [7^ 7l) 1 + Tr (7^^ [7^ 7l 7'^) 7. 
+ Tr (7^ [7^ Y] (t'^') ^.x 

+ Tr (7^^ [7^ 7l 7'^7') 7'7. + Tr (7^ b\ Y] Y) 7'} • 

Using the trace identities for gamma matrices (B.34), we note that traces of 
an odd number of gamma matrices vanish. For the two nonvanishing traces 
we obtain 

Tr (7^^ [7^ 7"] 7^^) = 4 {g^'g'" - g^'g"" + g^'^g'"' - g^'g'"' 

^g^'g'" -g^^g"') 

= 8 g^^'g"" - 8 gf'Pg''^ 

and 

Tr (7^^ [Y, Y] YY) = - 4 ? e^^'f" + 4 i e^^"" = -8 i e^""". 
Thus, we finally obtain the relation 

Y [Y, Y] = 2 g^'^Y - 2 g^^^Y - 2 i e^'^P'^Yla ■ (B.40) 
Using the algebra (B.9), we can additionally derive the following relation: 

[Y, Y] 7^ = 2 g^^Y - 2 g'^'^Y - 2 i e^^'^'^Yla ■ (B.41) 
From both of these relations we can easily derive the following identity: 

^^^/^-p- = - zY {YYY - g^^Y + g'^'Y - g^'Y) ■ (b.42) 

Gravitino Polarization Tensor Due to the Lorentz structure of the spin 
sum for the gravitino, it must be a linear combination of the following ten 
tensors [72]: 

PuVv 

m^/2gi,u, m^/2lt>,lv, 7mP^, IvPf,, 
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where we have added factors of 7713/2 to get the correct dimension of the spin 
sum according to equation (2.50). Therefore we have 



P^j,v{p) = A 1713/2 g^u + B 1713/2^ f,-f^ + C '^i + D '^i^p^ + E 



m3/2 



"^3/2 "^3/2 "^3/2 

Analogous to the case of Dirac spinors, we get two polarization sums for the 
positive and negative frequency solutions. The coefficients can be uniquely 
determined using the constraints from the Rarita-Schwinger equations for 
the mode functions (2.44), (2.45) and from the projector property of the 
polarization tensor (2.50). 

For the construction of the positive frequency spin sum, application of 
(2.44) gives 

A = F, B = G, C = H, D = I, E = J 
and application of (2.45) gives the further relations 

A + 4:B + D = 0, 4:C + 2A + E = 0, D - A-2B = 0, E - 2C = . 
The polarization tensor then becomes 

P;1'Ap) =A{p + m3/2) \gf.u-\l l^lu + - + 2 ^ 

I -J \ ^^3/2 ^^3/2 ^^3/2 

The projector property (2.50) finally gives A = —1. Therefore, the positive 
frequency polarization tensor of the gravitino reads 

P,l{p) = -{f + "^3/2) \ gt,v-\\ if.iu + 7a. - lu^^ + 2 

^ ^ [ 3 \^ m3/2 m3/2 

(B.43) 

which is equivalent to the form given in equation (2.51). 

Similarly, the negative frequency polarization tensor is determined. From 
(2.44), (2.45) and (2.50) we obtain 

A + F = 0, B + G = 0, G + H = 0, D + I = 0, E + J = , 
A + AB -D = 0, AG - 2A - E = , D + A + 2B = 0, E + 2G = 

and A = 1. Thus, the negative frequency polarization tensor of the gravitino 
reads 

P^^u{p) = -{f- ^3/2) \9^.u-^\ 7/.7^ + 7/.^^ - lu^^ + 2 

(B.44) 

which is equivalent to the form given in equation (2.52). 
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Appendix C 

Calculation of Decay Widths 



C.l Feynman Rules 

In this section we provide the necessary Feynman rules for the calculation of 
the gravitino decay widths. This set is taken from [29] and amended by the 
addition of the rules for the negative frequency solution of the gravitino. 

Gravitinos are depicted as double solid lines, chiral fermions are drawn 
as single solid lines and scalars are drawn as dashed lines. Gauge bosons 
are represented by wiggled lines and gauginos by solid lines with additional 
wiggled lines. 

The continuous fermion flow is independent of the fermion number flow 
carried by fermions and sfermions and also of the momentum direction. 

External Lines 

The momentum p flows from the left side to the right side for the external 
lines shown below. 
Scalar particles: 

= 1. 

« = u'{p) , 

- = u'{p), 

- = v'{p) , 

♦ = ^'{p) . 



Gauginos and matter fermions: 
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Gravitinos: 



Propagators 



The momentum p flows from tlie left side to the right side for the propagators 
shown below. 



Scalar particles: 

i 



— m\ 



5 



Matter fermions: 



z • > • J 



Gauginos: 



z {-f + m^) ^. . 
p2 — 



a i r\y^\y^\j t h = — ^ ^ d . 



p2 — 



Massless and massive gauge bosons: 



a, /i 



b,u = J- 



P 



' — m\ 

Gauge Vertices 

The momentum p flows into the vertex for the vertices shown below. 
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C.l. Feynman Rules 



Gravitino Vertices 

The momentum p flows into tlie vertex for tlie vertices sliown below. 



/ 



/ 



PI 



/ 



/ 



/i 



/i 



J 



a,P i 

/ 

/ 



[I 





J 



a,P i 

/ 

/ 




'PI 



/i J 
b,p,p 



/i 




y(«)p » p 



4M, 



PI 



4 Ml 



PI 



b, V 



P 



a 

c,p 



6, u 



717" 



P 



a 

c,p 




._|_y(o).V|7",7l 
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C.2 Kinematics 

Scalar Products 

Two-Body Decays Starting from 

P = Pl+P2, (C.l) 

we multiply by the three different four-momenta to obtain 

(P-P) = (P-Pl) + {P-P2), 

{P ■ Pi) = {Pi ■ Pi) + {Pi ■ P2) , (C.2) 

{P-P2) = {Pl ■P2) + {P2-P2) ■ 

Using [p ■ p) = and {pi ■ pi) = mf yields the relations 

{P ■ Pl) 



M2 


+ ml 


- ml 




2 






— m^ 


+ ml 




2 




M2 


2 

— m| 


-ml 



{P-P2) = - ^^1^, (C.3) 

{Pi ■ P2) 



Decay Widths 



The partial decay width for the decay of a particle of mass M into n particles 
is given by [17] 



dT=^^\Mfd^n{p;Pu...,Pn), (C.4) 

where d^n is an element of the n-body phase space, given by 



i=l 



A(27r)3 2E, 



Two-Body Decays In this case, the partial decay width becomes 
or if we average over the spin states 
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C.3. Gravitino Decay Channels 



The momentum of the final state particles is determined by four-momentum 
conservation: 

1^1 = ^^(M2 - (mi + ma)') (M^ - (mi - ma)') 



2M ^ ' ^ y M2 - (mi +m2)' 



C.3 Gravitino Decay Channels 

In this section we compute the dominant decay channels of the LSP gravitino 
in the framework of bilinear i?-parity breaking. The relative decay widths of 
the different decay channels determine the branching ratios for these chan- 
nels. Therefore, this is a crucial input for the prediction of spectra of gravitino 
decay products. 



At tree level the decay of the gravitino into the boson and the tau neutrino 
gets contributions from two Feynman diagrams: 




Apart from the diagram with the abelian 3-vertex and subsequent zino-tau 
neutrino mixing through the tau sneutrino VEV, we have a diagram with 
the non-abelian 4-vertex that reduces to an effective 3-vertex using the tau 
sneutrino VEV from bilinear i?-parity breaking. Analogous to the decay of 
the gravitino into the photon and the tau neutrino, the amplitude for this 
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process reads 



iM = — u 




-(g) (1 + 7^)7^ [k 7l^r(P)^i*(^) 



7^ 7iv^r(p)^i*(^) 



Similar to the procedure in Section 3.2, we rewrite some prefactors and in- 
troduce the zino-zino mixing parameter 



We note that the zino-zino mixing parameter is real- valued. Thus, the squared 
amplitude for this process becomes 



" " 256 Ml, ~~^J 

■ {Ul^ Tr [(^ + m.) (1 + 7^) Y [I Y] P^M [Y . P\ Y (l - 7^)] 
+ 2 mz U^^ Tr [(^ + m,) (l + 7') 7^ % Y] P^M YY (l - 7^)] 
+ 2 mz U^^ Tr [(^ + m,) (l + 7') 7^7^^^^+ (p) [7^ f\ Y (l - 7^)] 
-f 4m| Tr [(^ + m.) (l + 7^) YYP^lM YY (l - 7^)] } 

= - Y^ll^ {Uh Tr [i (1 + 7^) Y [t, Y] P^P) [7p, ^] Y] 

+ 4 mz Uzz Tr (1 + 7') Y Y\ P^M lpY\ 

-f 4 Tr (1 + 7^) Y (mlYP^P) 7p - ^ P^P) ^) 7^] } 



• { (2 (p ■ fc)^ {p-q) + 2 ml/2 {p ■k){k-q)- m|m^/2 (P " ?)) 
- 2 mzm3/2 (2 (p ■ /c) (p ■ g) + ml/^ {k ■ q)) U^z 




a=l 



(C.9) 
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In the last step we already replaced squared four-momenta by the correspond- 
ing squared particle masses. The other scalar products of the four-momenta 
in this process are given by 



(p-k) 

(p-q) 

(k-q) 



2 

"^3/2 


+ m| 


-ml 




2 




"^3/2 


-m| 


+ ml 


"^3/2 


2 

-m| 


-ml 


2 



Therefore, we obtain for the squared amplitude 

\M. I ~ — — J , I (3 mt/o — m\m\,o — mtml,o — m% + ml 

+ ("^3/2 -3 m|) +ml (3rr4 -5m^/2)) 
- 4 mzm3/2 (2 mj/^ - mlmj/^ - m% - ml 
+ ml {2ml -ml/2)) Uzz 

+ 2 ^3/2 + 9 m\m\/2 - 9 m%ml/2 - m% + ml 

- ml {ml/2 + 3 m|) + ml (3 m| + 10 m\ml/2 - ml/2) ) } . 

Discarding the terms proportional to the negligible neutrino mass and taking 
into account the correct phase space factor, the width for the decay of the 
gravitino into the Z° boson and the tau neutrino becomes 



(C.IO) 



where the kinematic factors /S^, fz, jz and hz are given by 



^ m 



2 ' 

"^3/2 



2 1 m% 



3 7713^2 3 m\/^ 



1 mi 
2^3/2 



(C.ll) 



mx 



/ix = 1 + 10 ^ + 



2 4 
^3/2 ^3/2 
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The conjugate process il)z/2 Z^J^t, that produces tau antineutrinos, has the 
same decay width: 




(C.12) 

V'3/2 W+T- 

Two tree-level Feynman diagrams contribute to the decay of the gravitino 
into the boson and the r~ lepton: 




Using the rotation of the charged gauge eigenstates into the charginos (cf. 
Section 2.1) 

2 2 

^^ = Y.^W^X and W- = Y,U^-J~, (C.13) 

and the relation between the gauge eigenstates and mass eigenstates of the 
electroweak gauge bosons (2.8) together with the generators of the corre- 
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spending gauge groups (2.58), we can write for the amplitude of this process 
^M = - r (g) ^V2g (z>.) - ^ ^) E 



q^-mi± y 4Mpi 



8Mpil^2^ 

Rewriting some coefficients and introducing the wino-wino mixing parameter 

a=l A-a 

that is real-valued, the squared amplitude for this process becomes 
,2 C f 



l-^l = " 128M|; V^'"^ ml, 

■ {U^^ Tr [(^ + m^) (1 + 7^) Y [I Y] P^.ip) [7^ ^] Y (l - 7^)] 
+ 2 m^. f/,^H^ Tr [(^ + m.) (l + 7') 7^ 7I (p) 7'^7^ (l - 7^)] 
+ 2 mw U^^ Tr [(^ + m.) (l + 7') 7^^+ (p) [7^ ^] 7^^ (l - 7^)] 
+4ml, Tr [(^ + m,) (l + 7^) 7^^+ (p) 7^ (l - 7^)] } • 

The traces are completely the same as in the previously discussed decay. Only 
the prefactors and the kinematics differ. Adopting the correct prefactors, 
neglecting rrir and exchanging mz with mw in the kinematic factors, the 
width for the decay of the gravitino into the boson and the r~ lepton 
reads 

(C.15) 

where the kinematic factors Pw, fw: jw and hw are the same as those defined 
in equation (C.ll). 
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The conjugate process 1P3/2 ^ W r+ has the same decay width 



w- w- 




(C.16) 



h Vt 



The decay of the gravitino into the hghtest Higgs boson and the tau neutrino 
also gets contributions from two tree-level Feynman diagrams 



h 





In the MSSM decoupling limit the mixing angle a becomes equal to /3 — 7r/2 
(cf. Section 2.1). Therefore we have 



cos a 
sin a 



sin/3 , 
— cos /3 



(C.17) 



and the decomposition of the gauge eigenstates into the mass eigenstates for 
the neutral Higgs bosons (2.11) becomes 



1 /sin /3 — cos /3 



ycos/3 sin/3 
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(C.18) 



C.3. Gravitino Decay Channels 



The inversion of this equation gives for the hghtest Higgs boson 

h = V2{Re [H^] sin /3 + Re [H^] cos /3) . (C. 19) 

Apphcation of this relation to the Higgs-slepton mixing terms in equation 
(3.4) together with the relation for the sneutrino VEV (3.6) results in the 
coupling 

= -^ml ^ulh + h.c. (C.20) 

between the tau sneutrino and the lightest Higgs boson, that appears in the 
first Feynman diagram. For the second diagram we use the vertices of the 
gravitino with the Higgs gauge eigenstates Re[if° and the corresponding 
higgsino states ^. Then we use the rotation of the Higgs gauge eigenstates 
Re ^] into the mass eigenstates h and H, which is given above in equation 
(C.18). Since we assume that the heavier Higgs bosons are kinematically not 
accessible, we only consider the final states with the lightest Higgs boson. 
Then we use, analogous to the other gravitino decay modes, the rotation of 
the neutral gauge eigenstates into the neutralinos (cf. Section 2.1) 

4 

K<i= E ^PLxl, (C.21) 

i, a=l 

and the relation between the gauge eigenstates and mass eigenstates of the 
electroweak gauge bosons (2.7) together with the generators of the corre- 
sponding gauge groups (2.58) to write for the amphtude of this process 

+ M''(g) iV2 [g cos dw - g'Y^, sin dw^ Pr 

* \i,j,a=l Xa / 

I COS/3 ( A Hl + rn^ 

\i,j,a=l Xa 



V2Mpi 
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igz (j>r) 



+ 



sm 



4t 



4 Q* Q 



.a=l 



4a/2Mpi 

Introducing the higgsino-zino mixing parameters 

4 C* C - 



4: Q* Q 

+ COS p > 



a=l 



^0 z - mz 

u, a 



E 



Q!=l 



(C.22) 



and rewriting some prefactors, we obtain for the squared amphtude 



\M\ 



64M2j 



2 



+ sin /3 f/ fro V + cos /3 f/fi 



IV [(^ + m.) (1 + 7^) P+ (p) ^ t'^ (1 - 7') ] 



32M2j 



+ sin PUfjoz + cos pUfjoz 



+ sin /3 [7^02 + cos pUjjoz 



■ {{P ■ kf (p-q)- m\m\i^ (p • ?)} ■ 

Here, we aheady substituted the squared four-momenta by the corresponding 
squared particle masses. The other scalar products of the four-momenta in 
this process are given by 



Therefore, we get for the squared amplitude 



2 

"^3/2 




-ml 




2 




"^3/2 


-ml 


+ ml 


2 

"^3/2 


2 

-ml 


-ml 


2 



M 



24 m^/^Mli 



m7. 



m-, — mi 



-\-sm^Ufjoz + cos PUfjoz 



r/2 2\3 2/4io22 o4\ 

I ("^3/2 -mh) -m^ (m3/2 + 2 m^mg/s - 3 m^j 
- mt {ml/2 + 3ml) +miy 
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Discarding the terms proportional to the neghgible neutrino mass and taking 
into account the correct phase space factor, the width for the decay of the 
gravitino into the hghtest Higgs boson and the tau neutrino becomes 



C ^3/2 

3847r 



f3t 



mf, — m 



+ sin (3 Ujjo z + cos (3 Uj^o ^ 



h 



(C.23) 

where the kinematic factor is given in equation (C.ll). 

The conjugate process 1^3/2 hur, that produces tau antineutrinos, has 
the same decay width: 





r (V'3/2 h Vr) 



e 



2 772^ 



3847r M, 



3/2 oA 
T Ph 



PI 



2 

mf, 



m~ — mt 



+ sin (3 Ujjo z + cos (3 Uj^o^ 



(C.24) 
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Appendix D 
PYTHIA Simulation 



In this appendix, we present an example of the PYTHIA simulation used to 
obtain the spectra of the fragmentation products. 

PYTHIA is a Monte Carlo event generator for collider simulations. There- 
fore we have to mimic a gravitino decay in a collider setup. To implement 
the gravitino in the simulation, we use the optional fourth fermion generation 
and redefine the properties of the particles. 

The process starts with the exclusive production of a Z' boson with mass 
mzi — 2 777.3/2 in an e" collision with a center of mass energy slightly above 
twice the gravitino mass. We redefine the decay channels of the Z' boson, so 
that only the decay into v'^ v'^ is allowed. The fourth generation neutrinos v'^ 
and v'^ are chosen to represent gravitinos with mass 7723/2. In order to have a 
nonvanishing phase space for this decay, we chose above a slightly increased 
center of mass energy for the initial collision. 

Then we define a single decay channel for the v'^: Either into Zv^-^ W^r^ 
or hi^r- The introduction of only one decay channel at a time allows us to 
determine the spectra independent of the branching ratios for the different 
decay channels. Anyway, these can easily be added afterwards. 

We are interested in the stable final state particles, so we have to allow 
all decays of long-lived particles like muons, pions, kaons and neutrons. By 
default, PYTHIA regards these particles as stable, since they do not decay 
on collider timescales. 

In the end we print out a list with the energies of the final state neutrinos 
(only one flavor at a time). 

The used PYTHIA program is written in FORTRAN 77 and is compiled 
via the command 

g77 ' program '.f pythia-6.4. 15.f . 

This produces an executable file a. out. The event generation is started via 
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a . out > ' output ' . dat 



and the output data is stored into the file ' output ' . dat. The following code 
shows the program for the simulation of the decay and fragmentation of 
ip3/2 — > Z'^i't- and the corresponding conjugate process. 

C. . .Double precision and integer declarations. 
IMPLICIT DOUBLE PRECISION (A-H,0-Z) 
IMPLICIT INTEGER (I-N) 
INTEGER PYK,PYCHGE,PYCOMP 

C... EXTERNAL statement links PYDATA on most platforms. 
EXTERNAL PYDATA 

C. . .Common blocks. 
C. . .The event record. 

C0MM0N/PYJETS/N,NPAD,K(4000,5) ,P(4000,5) ,V(4000,5) 
C. . .Parameters . 

C0MM0N/PYDAT1/MSTU(200) ,PARU(200) ,MSTJ(200) ,PARJ(200) 
C. . .Particle properties + some flavor parameters. 

C0MM0N/PYDAT2/KCHG(500,4) ,PMAS(500,4) , PARE (2000) ,VCKM(4,4) 
C. . .Decay information. 

C0MM0N/PYDAT3/MDCY(500,3) ,MDME(8000,2) ,BRAT(8000) ,KFDP(8000,5) 
C ... Selection of hard scattering subprocesses . 

C0MM0N/PYSUBS/MSEL,MSELPD,MSUB(500) ,KFIN(2,-40:40) ,CKIN(200) 
C. . .Parameters . 

C0MM0N/PYPARS/MSTP(200) ,PARP(200) ,MSTI(200) ,PARI(200) 
C. . . Supersymmetry parameters. 

COMMON/PYMSSM/IMSS (0 : 99) , RMSS (0 : 99) 

C. . .Set gravitino mass. 
GMASS=150 . ODO 

C. . .Set number of events and CMS energy. 

NEVENTS=200000 

ECM=2*GMASS+0 . 003 
C...Note: ECM is slightly increased so that PYTHIA doesn't complain. 

C. . .Allow fourth fermion generation. 
MSTP(1)=4 

C. . .Set Z'O (32) mass to 2*gravitino mass. 
PMAS(32,1)=2*GMASS 
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C...Set nu'_tau (18) mass to gravitino mass. 
PMAS(18,1)=GMASS 

C. . .Allow only Z'O production e+e- — > Z'O. 
MSEL=21 

C. . .Allow decays of nu'_taus and Z'Os. 

MDCY(18,1)=1 
MDCY(32,1)=1 

C. . .Turn off all decay channels of the Z'O. 

END=MDCY (32,2) +MDCY ( 32 , 3 ) - 1 
C.Note: MDCY(32,2)=289, MDCY(32,2)+MDCY(32,3)-1=310 

DO 300 IDC=MDCY(32,2) ,END 

MDME(IDC,1)=-1 

BRAT (IDC) =0 
300 CONTINUE 

C. . .Turn on Z'O — > nu'_tau nu'_taubar decay channel. 

MDME(MDCY(32,2)+15,1)=1 

BRAT(MDCY(32,2)+15)=1 
C.Note: MDCY(32,2)+15=304 

C. . .Define decay channel nu'_tau — > ZO (23) and nu_tau (16). 
KFDP (MDCY (18 , 2) , 1) =23 
KFDP(MDCY(18,2) ,2) =16 

C. . .Set branching ratio nu'_tau — > ZO nu_tau = 1. 
MDME(MDCY(18,2) ,2)=0 
BRAT (MDCY (18, 2) )=1 

C. . .Turn off other decay channels of the nu'_tau. 

END2=MDCY (18,2) +MDCY ( 18 , 3) - 1 
C.Note: MDCY (18, 2) =150, MDCY(18,2)+MDCY(18,3)-1=152. 

DO 200 IDC=MDCY(18,2)+1,END2 

MDME(IDC,1)=-1 

BRAT (IDC) =0 

MDME(IDC,2)=102 
200 CONTINUE 

C. . .Add new decay channel for the neutron. 
MDCY(PYC0MP(2112) , 1)=1 
MDCY(PYC0MP(2112) ,2)=7000 
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MDCY(PYC0MP(2112) ,3)=1 



C. . .Define neutron decay channel n — > p e- nu_e. 
KFDP(MDCY(PYC0MP(2112) ,2) ,1)=2212 
KFDP(MDCY(PYC0MP(2112) ,2) ,2)=11 
KFDP(MDCY(PYC0MP(2112) ,2) ,3) =12 
MDME(7000,1)=1 
BRAT (7000) =1 

C... Increase neutron mass slightly so that PYTHIA doesn't complain. 
PMAS(PYC0MP(2112) ,1)=0.942 

C. . .Turn off initial state radiation. 
MSUB(1)=1 
MSTP(11)=0 

C. . .Allow (21) and force (22) all unstable particles to decay. 
MSTJ(21)=1 
MSTJ(22)=1 

C... Allow decays of muons(13), pions(211), K_L0s(130) and K+-s(321). 
MDCY(13,1)=1 
MDCY(PYC0MP(130) ,1)=1 
MDCY(PYC0MP(211) ,1)=1 
MDCY(PYC0MP(321) ,1)=1 

C. . .Allow neutron decay. 

MDCY(PYC0MP(2112) ,1)=1 

C. . .Uncomment to show a list of decay properties. 
C CALL PYSTAT(2) 

C. . .Uncomment to show a list of decay channels. 
C CALL PYLIST(12) 

C. .. Initialization of the event generation procedure. 
CALL PYINIT ( ' CMS ' , ' e+ ' , ' e- ' , ECM) 

C. . .Generation of events. 

DO 100 NEV=1,NEVENTS 
CALL PYEVNT 

C. . .Uncomment to show a summary of the event listing. 
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C CALL PYLIST(l) 

C. . .Print the energies of final state neutrinos and antineutrinos : 
nu_e=12,-12, nu_mu=14,-14, nu_tau=16,-16. 
DO 400 1=1, N 

1F(K(I,1) .EQ.l.AND. ((K(I,2) .EQ.16) .OR. (K(I,2) .EQ.-16))) 
THEN PRINT *,P(1,4) 
END IF 
400 CONTINUE 
100 CONTINUE 

END 
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